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Abstract Freeze-out of the gas phase elements onto cold grains in dense interstellar and circnm- 
stellar media builds up ice mantles consisting of molecules that are mostly formed in situ (H 2 O, 
NH 3 , CO 2 , CO, CH 3 OH, and more). This review summarizes the detected infrared spectroscopic 
ice features and compares the abundances across Galactic, extragalactic, and solar system envi¬ 
ronments. A tremendous amount of information is contained in the ice band profiles. Laboratory 
experiments play a critical role in the analysis of the observations. Strong evidence is found for dis¬ 
tinct ice formation stages, separated by CO freeze out at high densities. The ice bands have proven 
to be excellent probes of the thermal history of their environment. The evidence for the long-held 
idea that processing of ices by energetic photons and cosmic rays produces complex molecules is 
weak. Recent state of the art observations show promise for much progress in this area with planned 
infrared facilities. 
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1 PERSPECTIVE 


(Interstellar) ice: 

solid state pure 
volatiles (H 2 O, CO, 
etc.) or mixtures 
thereof, with, per¬ 
haps, refractory 
inclusions, located 
on interstellar dust 
grains. 


The condensation of gas phase elements o n cold dus t parti cles has long been consid¬ 
ered an important process in the Universe. iLindbladl (jl935l i studied the condensation 
of metals to form meteoritic materi al as the starting poi nt of planet formation. Con¬ 


sidering the elemental abundances, van de Hulst ( 1946l i realized that the same pro¬ 


cess must lead to the formation of H 2 O, CH 4 , and NH 3 ices on the grains. Chemistry 
on cold grain surfaces has since emerged as a vital route to formation of new molecules 
in dense clouds and circumstellar en velopes and disks, including species that cannot 
readi ly be formed in the gas phase (jTielens &: Hagenlll982l : iHerbst &: van Dishoeck 


2 OO 9 I and references therein); surface chemistry thus controls the composition not 


only of the icy mantles but also of the molecular gas desorbed from the grains in 
the vicinity of cloud edges, shocks, and the radiation fields of Young Stellar Objects 
(YSOs). 

Other chemical pathways to new interstellar ice species are observationally not 
as well established. Laboratory studies have shown that ice heating creates new 
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molecules by purely thermal reactions ( Schutte et ah . 19931 : Theule et ah . 2013l l or 
through the increased mobility of radicals created by ultraviolet (UV) photons and 
cosmic ray (CR) particles. The latter has been a central topic in laboratory studies 
as a means to form complex molecules from simple ices. Heavy processing of ices 
forms organic residues (sometimes referred to as “yellow” and “brown s tuff”) that 
resemble the insoluble orga nic component of carbonaceous meteorites ([Greenberg 
1973 : Greenberg et al. 19951 . and references therein), and the products may incl ude 


biologically relevant species fe.g.. iBernstein et al.lll995l : iMuhoz Garo et al.ll2004^ . 

Ices are an important reservoir of the elements in dense (> 10^ cm“^) environ¬ 
ments; up to 60% of the oxygen not included in silicates is incorporated in known 
ice species. The processes that are key in their origin and evolution are thus di¬ 
rectly relevant to key questions on the organic and volatile reservoir of regions of 
star and planet formation. What is the origin of the volatiles and organics of the 
solar nebula, possibly delivered to the early Earth as the building blocks of life? 
There is no clear consensus on whether these materia ls survived the voyage from 
circumstellar envelo 


2014 : ICleeves et al 


res in to protoplanetary disks (e.g.. IVisser et al. 


2009; 


Brownlee 


2014 ). Gomets, believed to contain the most pristine remnant 


ices from the formation of the solar system le.g.. lMumma fc Charnlevll201ll ). are key 
in answering this question. 

In the broader astrophysical context, ices pla y important roles a s well. For exam¬ 
ple, they speed up the grain coagulation process ( Ormel et al. . 20091 ). Also, questions 
on the importance of ice processing are directly related to questions about the origin, 
composition, and physical properties of interstellar dust. Energetically processed ices 
might be a sou rce of refractory interstellar dust, in particular in star forming regions 
(jLi et al.ll20l4 and references therein) as opposed to dust formed in the envelopes of 
evolved stars. 

The initial searches fo r ices targeted t h e 3.0 pm HpO band in line s of sight dom¬ 


i nated by diffuse dust (jPanielson et al.l . Il965l : iKnacke et al 


196 


3 ). 


Gillett et al 


( I 975 I ) noted that the lack of ices in these diffuse clouds could not be explained by 
elevated dust temperatu res alone, revealing the im portance of the non-thermal pro¬ 
cess of photodesorption (jWatson &: Saloeterl. Il972l). The 3.0 um f eature was finally 
detected toward the Orion BN/KL region ( Gillett &: Forrest . 19731 ). Its position and 
width were in reasonable agreement with a model of small, pur e H 91 O ice spheres. 


but significant deviations were observed as well ( Merrill et ah . 19761 ): short and 


long-wavelength wings were attributed to scatter i ng by large gra i ns and interactions 


with NH 3 mixed in the ice (jMerrill et al.l . 119761 : iKnacke et al.l . Il982l : ISmith et al. 


198 


^). The latter sup ported the concept of ice mixtures (sometimes called “dirty 


ice”) first proposed by van de Hulst ( 19461 ). Most ice features discovered thereafter 


show stunningly complex profiles as a result of the ice composition, structure (ther¬ 
mal history), and grain shape and size. For several bands (those of H 2 O, CO, and 
CO 2 in particular) this is well understood, and their diagnostic value of the physical 
conditions and the nature of astronomical sources is a main topic of this review. 

This is a review of observational studies of interstellar and circumstellar ices. 
The composition of the ices and evidence for formation and thermal and energetic 
modification processes are reviewed across Galactic and extragala ctic environments . 
Previous reviews focused on rnassive YSOs (MYSOs hereafter; iGibb et al.l 12004 : 
Boogert fc Ehrenfreund 2004b : Dartois 20051 ) and on low mass YSO s (LYSOs here¬ 


after) and background stars tracing nearby clouds ([Oberg et al.l . l201ll ). In particular. 


Energetic pro¬ 
cessing of ices: 

exposure to ionizing 
radiation, breaking 
bonds, and creating 
radicals and ions. 
Subsequent chemical 
reactions may lead 
to more complex 
molecules. 

Photodesorption: 

non-thermal desorp¬ 
tion of a particle 
from a grain sur¬ 
face after a photon 
hit to itself or an 
underlying particle. 


Dirty ice: ill- 

defined term some¬ 
times used to 
indicate that in¬ 
terstellar ices are 
mixtures, contain 
processed residues, 
or are present on 
grains. 
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the study by Oberg et al. is expanded here to include the Galactic Center (GC) re¬ 
gion and extragalactic environments, as well as to a full inventory of observed features 
and identifications. Solar system ices are not reviewed in detail, but the interstellar 
medium (ISM) and circumstellar ice prope rties are compared with those of comets 
as reviewed by Mumma Sz Charnlev ( 201lh to address the question of the origin of 
solar system volatiles and organics. 


2 ICES ACROSS THE ELECTROMAGNETIC SPECTRUM 


Ices in interstellar environments are traced almost exclusively by their molecular 
vibrational transitions in the near- to far-infrared. In general, the 1-3 pm wave¬ 
length region contains the combination and overtone modes, the 3-6 pm region the 
stretch vibrations, the 6-30 pm region the bend and libration vibrations, and the 
25-300 pm region the torsional and inter-molecular (lattice) modes. Some ices also 
have strong electronic tra nsitions in the UV, e.g. in the case of H 2 O near 0.14 pm 
(jWarren &: BrandtI . l2008l ). These UV features have not yet been observed in the 
ISM, due at least in part to very high UV extinction in candidate lines of sight. 
They may be detectable in absorption near cloud edges (at low Ay values; iGoebel 
1983l l or through scattering. Finally, simple ices are transparent outside their vibra- 
tional modes, but the continuum opaci ties increase significantly for processed ices 
( Jenniskens . 199^Brunetto et ah . 2006l l and ices with amorphous carbon inclusions 


(jPreibisch et al 


199. 


W- 


The vast majority of interstellar and circumstellar ice studies target the stretch 
and bend mode vibrations in the 3-16 pm range, because below 3 pm the signals are 
weak as a result of dust continuum extinction (in contrast to solar system studies 
which often observe the combination modes at 1-2 pm), and above ~30 pm the 
availability and capability of instrumentation are limited. Figures [T] and [5] show 
a selection of ice features detected outside the solar system and Table [U gives a 
comprehensive summary of the typical peak positions, widths, and identifications. 
Nearly all features are securely detected, but not all are positively identified, as 
indicated with questions marks and further discussed in 116.11 The widths vary 
dramatically, from 1.5 cm“^ (0.0034 pm) for ^^GO to 400 cm“^ (0.38 pm) for the 
H 2 O ice bands. Also, many band profiles (e.g., those of GO and GO 2 ) contain 
substructures at the level of a few cm“^. Thus, ice studies benefit significantly from 
the availability of medium resolution infrared spectrometers {R = A/AA >500). For 
many features, the positions and widths vary considerably across objects. These 
variations are often well studied, and have led to empirical decompositions. The 
components are listed, and are further discussed in @ 

Indirect constraints on the ice composition are provided by observations of gas- 
phase rotational lines at (sub)-millimeter wavelengths. First, low gas phase abun- 
dances (depletions) of s pecific molecules in clouds imply their presence in the ices 
(jBergin &: Tafallal . 120071 1. Second, gas phase species observed in hot core regions 
surrounding YSOs originate from sub l imated ices, either directly or after gas phase 

Thir d, some 
). These 


chemical reactions ( Blake et al.l . 119871 : iHerbst &: van Dishoeck . 2003 1 


2009| l 

3 t al.l . 


gas near cloud edges originates from photodesorbed ices (jOberg et al 
results are not reviewed here, but are quoted when relevant. 
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Figure 1: Overview o f the strongest ice and dust features in the MYSO AFGL 7009S 

The calculated spectrum of pure H 2 O ice spheres at 10 K is 


(jPartois et al. 


shown (dashed line) to indicate the multiple H 2 O bands. 


2.1 Ice Mapping 

Most ice features are detected as pure absorption bands against infrared continuum 
point sources. Studies of the spatial variations of the ice properties are therefore 
relatively rare, yet very powerful. They often rely on the presen ce of many point 


sourc es (usually infrared-bright giants) behind clouds (e.g., Taurus: iMurakawa et al 


200Cli. The envelope of a C lass 0 YSO was mapped using background Class II LYSOs 
(IPontoDPidan et al.l . l2004l b and of a dense core using Class I/II LYSOs (|Pontoppidanl . 
200y). Ice mapping is however also possible if the background emission is extended 
by scattering frorn dust in any disk, e nvelo pe, or outflow cone (jHarker et al.l . 119971 : 
Spoon et ah . 2003 : Schegerer &: Wolf . 201 oI l In addition, scatterin g by large grains 


enha nces the short-wavelength wing of the 3.0 pm H 2 O band f ll4.lt iPendleton et al 


199oI l which enables mapping of ices in reflection nebulae. This has rarely been 
done. Similarly, the lattice modes (~ 25-300 pm; 114.51) are suited for ice mapping. 
They are the only ice bands that may appear in emission, because, in contrast to 
the intra-molecular modes, their excitation energy is below the binding energy of the 


ices. 























































6 


Boogert, Gerakines, Sz Whittet 



Figure 2: Compilation of weak ice absorption features in the 3.3-4.0 (top), 7.0-7.9 
(middle) and 8.0-10.4 (bottom) gm regions. Identifications are labeled, and tentative 
identifications are fo llowed by questio n marks. Data for A FGL 2136 were taken from 


Dartois et al. ( 

2003 

, for W 33A from 

Schutte et a 

1 . (lOOflI'l. and for IR AS 15398-3359 

from 

Booeert et al. 

(2008ll and 

Bottinelli et al. (2 

0101 . 














































Table 1; Ice Features Detected in Interstellar and Circum- 
stellar Environments 


^center 

fim 

^center 

cm“^ 

FWHM 

cm“^ 

Identification 

Reference®" 

2.27 

4405 

60 

CH 3 OH combination 

1 

2.697 

3708 

4 

CO 2 combination 

2 

2.778 

3600 

5 

CO 2 combination 

2 

2.96 

3378 

~ 40 

NH 3 stretch? 

3,4 

3.0 

3333 

~ 390 

H 2 O stretch (multiple components)*’ 


—2.90 

3448 

~ 120 

— scattering 

5, 6 

—3.05 

3278 

~ 315 

— amorphous small grains 

7, 8 

—3.10 

3225 

~ 220 

— crystalline small grains 

8 

-3.3 

3030 

~ 270 

— scattering and NH 3 .H 2 O hydrate 

9, 10, 16 

3.25 

3076 

70 

PAH C-H stretch? NHJ combination? 

11, 12, 13 

3.325 

3008 

15 

CH 4 stretch 

14 

3.47 

2881 

80 

NH 3 .H 2 O hydrate 

15, 4 

3.53 

2832 

30 

CH 3 OH C-H stretch 

17 

3.85 

2597 

~ 40 

CH 3 OH combination 

15 

3.94 

2538 

~ 40 

CH 3 OH combination 

18, 15 

4.07^ 

2457 

~ 120 

HDO 0-D stretch amorphous 

19 

4.13'= 

2421 

~ 60 

HDO 0-D stretch crystalline 

19 

4.27 

2341 

~ 20 

CO 2 stretch 

20 

4.39 

2278 

9 

*^C 02 stretch (multiple components) 


—4.382 

2282 

3 

— pure CO 2 

21 

—4.390 

2278 

8 

— polar, H 2 O/CO 2 >1 

21 

4.5 

2222 

270 

H 2 O combination 

22 

4.62 

2164 

29 

“XCN” (multiple components) 

23 
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'^center 

^center 

cm“^ 

FWHM 

cm“^ 

Identification 

Reference®' 

—4.598 

2175 

15 

—OCN“ C=N stretch (apolar)? 

24, 

25 

—4.617 

2166 

26 

—OCN“ C=N stretch (polar) 

24, 

25 

4.67 

2141.3 

3-9 

CO stretch (multiple components) 

26 


—4.665 

2143.7 

3.0 

— apolar, C 02 /CO> 1 ? 

44, 

24 

—4.673 

2139.9 

3.5 

— apolar, > 90% CO 

23, 

24 

—4.681 

2136.5 

10.6 

— polar, CH 30 H/C 0 > 1 

23, 

24, 27 

4.779 

2092.4 

1.50 

^^CO stretch (apolar) 

28 


4.90 

2040 

22 

OCS C-0 stretch 

18, 

29 

5.83 

1715 

35 

H 2 CO C-0 stretch 

30 


6.0 

1666 

160 

multiple components 

31 


— 6.0 

1666 

~ 130 

—H 2 O bend 

31 


—5.84 

1712 

72 

—HCOOH C-0 stretch? 

30, 

31, 32 

—6.18 

1618 

111 

—NH 3 bend, H 20 /C 02 = 2 , PAH C-C stretch 

30, 

32 

6.0 

~ 1666 

~ 330 

organic residue or an-ions? 

33, 

32 

6.85 

1459 

85 

multiple components 

31 


—6.755 

1480 

43 

—CH 3 OH C-H deformation, NH^ bend? 

30, 

32 

—6.943 

1440 

61 

—NH| bend? 

30, 

32 

7.24 

1381 

19 

HCOOH, CH 3 CH 2 OH, C-H deformation? 

34 


7.41 

1349 

15 

HCOO“, CH 3 CHO, C-H deformation? 

34 


7.63 

1310 

28 

SO 2 stretch? 

35 


7.674 

1303 

11 

CH 4 deformation 

36, 

37 

8.865 

1128 

15 

CH 3 OH C-H 3 rock 

38 


9.01 

1109 

50 

NH 3 umbrella 

39 


9.74 

1026 

30 

CH 3 OH C-0 stretch 

38 


13.6 

735 

~ 220 

H 2 O libration 

40 
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-^center 

fim 

^center 

cm“^ 

FWHM 

cm“^ 

Identification 

Reference*^ 

15.1-15.3 

653-662 

16-28 

CO 2 bend (multiple components) 

20 

—15.31 

653 

22 

— polar, C 02 /H 20 < 1 

24 

—15.10 

662.2 

11 

— apolar, C 02 /CO~l 

24 

—15.15 

660.0 

2.5 

— dilute, C 02 /CO<< 1 

24 

—15.11-bl5.26 

655.3-b661 

19 

— pure CO 2 

24 

—15.4 

648 

10 

— shoulder, CH 3 OH/CO 2 >> 1 

24 

44 

227 

~ 100 

H 2 O lattice (amorphous) 

41, 42, 43 

44 

227 

~ 50 

H 2 O lattice (crystalline) 

41, 42, 43 

62 

161 

~ 50 

H 2 O lattice (crystalline) 

41, 42, 43 


^ reference to a n 

ublication contributine 

stronelv to its identification or profile characterization: 1 : 

Taban et ah 

2003], 

2: I 

Ceane et al. 2 

001 a 

3 

Ghiar et al. 

2000, 4: 

Dartois & d 

Hendecourt 

200 

1, 5: 

Knacke et al. 

198^ 

!, 6 : 

Pendleton et al. 

1990. 7: Gillett & Forres 

t 1973. 8 : 

Bmith et al. 

1989, 9: 

Merrill et al. 

1976, 10: 

inacke & McGorklel ! 

-987, 11: 

Selleren et al. 

I994I. 12 : 

Hardeeree-Ullman et al. 

2014 

, 13: S 

chutte & Khanna 2003 

14: 

Booeert et al. 

2004a. 15: 

A-llamandola et al. 


1992 

16: 

20 : 

, 24: 

Dartois et al. 

200: 

2. 17: Grim et a 

. 1991, 18: 

Geballe et al. If 

85, 19: 

A-ikawa et al. 

i 

2012 

Gerakines et 

R. 1999. 21: I 

looeert et al. 

2000 a, 22 : 

Booeert et al. 

200c 

)b, 2 c 

■. 

27 

Lacv et al. 

1984 

PontoDoidan et al. 

2003 

25: 

van Broekhuizen et al 

2005. 2f 

: Soifer et al. 

1979, 

GuDoen et al. 

2011 , 

28: 

Booeert et al. 

2002 s 

, 29: 

Palumbo et al. 

1997, 

30: 

Keane et al. 

2001b. 31: 

Schutte et al. 

199f 

, 32: 

Booeert et al. 

2008. 33: 

Gibb & Whittet 

2002. 34: S 

Ichutte et al. 

1999. 35: 

Booeert et al. 

1997. 36: 

Lacv et al. 1991. 37: Booffert et al. 1^ 

)96, 3 

8 : 

43 

Skinner et al. i 

992. 39: 

Lacv et al. 

1998 

. 40: Gox 

1989. 41: 

Omont et al. 

1990. 42: I 

vlalfait et al. 

1999, 

: Dartois et al. 

1998. 44 

: Booeert et al. 

2002 hl 


Contributions from the 0-H stretch mode in CH 3 OH and perhaps HCOOH and other species 
are expected, but are not detected as separate features. 

Tentative detection. 


CO 
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3 LABORATORY EXPERIMENTS 


Laboratory experiments have been crucial in the identification of the observed inter¬ 
stellar ice features (Tabled]), in understanding their profiles (SI), and in determining 
the chemical origin and fate of their carriers (SI)- An overview of laboratory work 
is beyond the scope of this review, but key aspects will be emphasized here and 
mentioned throughout the remaining sections. 

The peak positions, widths, and shapes of the observed ice bands are influenced by 
the composition and structure of the ices (crystallinity, density, and homogeneity) 
as a result of the attractive and repulsive dipole interactions between neighboring 
molecules. They are also influenced by the properties of the dust grains (size, shape, 
and the relative mantle and core volumes) and these effects can be calculated from 
the optical properties of the bulk materials. All of these influences must be taken 
into account in the proper assignment of carriers to the observed ice features and 
special care must be used in deriving physical and chemical information from the 
band profiles such as overall ice composition or temperature. 

The absorption band profiles of ices and ice mixtures are routinely measured 
by laboratory transmission and reflection spectrosco py at high vacuum (< 10 “^ 


mbar) under temperature-controlled conditions (e.g., Hudgins et al. 1993lf . 


These 

laboratory samples may also be exposed to UV photons or to particles, simulating 
molecule fo rmation and destructio n by energetic processes in interstellar environ¬ 
ments fe.g.. iGerakines et al.ll2001bl ). Optical constants are generally measured with 
a Kramers-Kronig analysis of a single spectrum using the ice’s known thickness and 
visible refractive index. F or CO and CO 2 , conflicting optical constants were re¬ 


ported (|Ehrenfreund et al.l . 1 199711 , and it was shown that the most accurate values 


are d erived using polarized ligh t at a range of sample thicknesses and incidence an¬ 


gles ( Baratta Sz Palumbd . 119981 1. Databases of laboratory transmission spectra and 
optical constants of interstellar ice analogs are listed in ’’Related Resources”. 

Infrared band strengths (^ can be determined directly from the ice’s complex re¬ 
fractive index, or from laboratory spectra of an ice of a known thickness that must be 
determined independently. The band strengths of a particular molecule vary with ice 
temperature and overall ice composition, although measurements under astrophysi- 
cal conditions are relatively scarce. The band strengths of CO and CO 2 were shown 
to vary by less than '-^ 20 % when measured at l ow temperatures, mixed together. 


or in an H 2 O ice matrix (jCerakines et al.l . Il995l l. but the strength of the 3.0 gm 


H 2 O stretch mode is reduced by an order o f mag nitude when H 2 O is highly diluted 


(van Thiel et al.. 

0 ^ 

T—1 

Ehrenfreund et al. 

ice bands are listed in, e.g.. 

Hudffins et al. 


1996 1. Int egrated strength s of co mmon 
( 1993 1 and Cerakines et al. ( 1995 1. 


Improvements in laboratory techniques have made it possible to measure fun¬ 
damental parameter s such as energy barr iers, reaction rates, and photodesorption 
yields. We refer to iLinnartz et al.l ( 201 il l for a recent review on this topic. Very 
briefly, this “bottom up” app roach h as been succes s ful in creating H 2 CO and CH 3 OH 
(jWatanabe fc Kouchil . I 2 OO 2 I I. CO 2 (|Ioppo1o et al.l . I 2 OIII I and other species on cold 
surfaces with quantitative results. These have been used in grain surface-chemistry 
models and have successfully reproduced observe d ice abundances under realistic 
astrophysical conditions (e.g., Cuppen et al. 201 ill . 
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4 ICE BAND PROFILES 


The profiles of the 3.0, 4.67, and 15.2 pm ice features have proven to be sensitive 
probes of the physical conditions and the nature of a wide range of astronomi¬ 
cal objects (^: dense clouds, Galactic and extragalactic YSOs, evolved stars, the 
Galactic Center region, and external galaxies. The lattice modes above 25 pm show 
great promise as probes of the ice composition and processing history. The 6 and 
6.8 ^m features, though not fully identified, were shown to have diagnostic value as 
well. All other bands are discussed in 3^1 


4.1 The 3.0 pm H 2 O Band: Crystallization, Large Grains, and Hy¬ 
drates 

The interstellar 3.0 ^m band (Fig. [3|) is primarily due to the 0-H stretching mode 
of bulk H 2 O ice, because of the high abundance of this species and because of 
the very high intrinsic band strength induced by the hydrogen bonding network 
( van Thiel et ah . 1957I T. The bands of H 2 O ice mixed with apolar species peak be¬ 
low 3.0 pm, and, lacking the hy drogen bonding network, their intrinsic strengths are 


an order of magnitude smaller (jEhrenfreund et all . 1 19^ 1 . These features of isolated 


H 2 O have not been detected in interstellar environments. 

The profile of the 3.0 pm band is not fully explained by absor ption by small pure 
H 2 O ice spheres (jMerrill et al.l . ll976l : lKnacke &: McCorkld . 119871 1. Excess absorption 
is always present between 3.3 and 3.7 pm (Eig. [3]), and sometimes near 2.9 pm. The 
strengths of the long wavelength wing an d the peak depth of the 3. 0 um band corre¬ 


late well. Some va riations were reported (jPontoppidan et al.l . l2004l : iThi et al.l . 120061 : 


Noble et al.l . 1201311 . suggesting a dependence on H 2 O abundance, but this needs fur- 
ther work. In any case, Mie scattering calculations show that large grains (~ 0.5 pm; 
Smith et al.lll989l : ^ might explain much but not all of the long wavelength wing. 
The remaind er of the absorption i s attributed to the 0-H s t retch ing mode in ammo¬ 


nia hydrates (jKnacke et al 


1982 


Dartois &: d’Hendecourtl. 


pm f eature is likely also a result of these hydrates f (16.1l 

200 ih . 


2001). The distinct 3.47 


Dartois fc d’Hendecourt 


The wing at 2.9 pm is not as commonly observed. Its strength increases with the 
observational beam size, capturing larger contributions of light sc attered off large 


ice-co ated grains in some objects (e.g., the Orion BN/KL nebula in iPendleton et al 


I 990 I ). The absence of this feature in a background star and i ts presence toward s 


YSOs and reflection nebulae further confirms this interpretation ( Smith et ah . 19891 ) 


The stretching mode of NH,-; contr i butes to this wing, but it is ra rely detected as a 


distinct feature (Ghiar et al 


200 c : Dartois fc d’Hendecourt . 200 ll) . 


The 3.0 pm band profile varies considerably among targets within and across 
different environment classes (Eig. [3|). The peak position and width toward the 
maiority of observed YSOs and all background stars are consistent with amorphous 
H 2 O ice ( Smith et ah . 19891 ). Narrower p rofiles that are shifted to longer wavelengths 


toward some MYSOs and a few LYSOs (jBoogert et al.l . l2008l ) reflect higher degrees 
of crystallinity attributed to thermal processi ng f()8.2l). An extr eme case is that 


of the edge-on low mass disk source d216-0939 (|Terada &: Tokunagal . l2012al ). which 
shows a peak at ~3.20 pm caused by large (~ 0.8 ^m) grains consisting of nearly 
pure crystalline H 2 O ice. 


Hydrogen bond¬ 
ing network: 

dipole attraction 
between the elec¬ 
tropositive H from 
one molecule with 
the electronegative 
atom (e.g., O or N) 
from another. 


Ammonia hy¬ 
drates: a hydrogen 
bonding network 
between H 2 O and 
NH 3 , often noted as 
H 2 O.NH 3 . 


Thermal process¬ 
ing of ices: in¬ 
crease in ice tem¬ 
perature, leading to 
sublimation, a rear¬ 
rangement of molec¬ 
ular bonds (crystal¬ 
lization and segrega¬ 
tion), or chemical re¬ 
actions. 
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Evolutionary Scenario 
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Figure 3: Observed 3.0 band profiles, highlighting the universally present H 2 O 
ice components and their variations. The left panel shows sightlines with, from top 
to bottom, increasing degrees of crystallinity as a result of heating above tempera¬ 
tures of ~ 70 K: the backgro und star 2MASSJ 21 240614+4958310 ( Boogert et ah . 
2011), the LYSO R NO 91 ( Boogert et ah . 20081 ). the Gal actic MYSO GL2136 


JPartois et ah . 2002), and the edge-on disk source d216- 0939 (Terada fc Tokunaga . 


OH 231.8+4.2 ( Dartois et al. 


2012al ). The right panel shows the SMG MYSO 03 (lOliveira et al 
ULIRG galaxy MIPS 2230 3 at z=2.28 


2013 ), the 


(jSaiina et al.l l2009l ) , and the evolved star 
200^ ) . For comparison, the smoothed spectrum of the 


background star is plotted over each spectrum (red). The background star itself is 
compared with pure, amorphous H 2 O ice particles (GDE shapes) at 10 K (green 
area). If instead a H 2 O:NH 3 = 10 :l spectrum is used, part of the long-wavelength 
wing is explained by NH 3 hydrates (NH 3 .H 2 O; red area). The black area is the 
residual due to large (~ 0.35 /um) grains. 


Polar ices: ices 

dominated by 

molecules with high 
dipole moments 
(e.g., H 2 O, GH 3 OH), 
which have relatively 
high sublimation 
temperatures. 


4.2 The 4.67 /im CO Band: the Most Volatile Ices 

The 4.67 gm band consists of three components as a result of solid CO in different 
molecular environments. Gaussian representations of these components fit all ob¬ 
served CO ice bands toward YSOs and background stars (Tabled! Fig. dj), and their 
relative strengths may b e used to derive th e ice c omposition and dust temperature 
along the line of sight ( Pontoppidan et al. . 20031 ). A narrow feature centered on 


4.673 /im and a much broader long wavelength wing peaking at 4.681 /im are tradi¬ 
tionally attributed to CO emb edded in low (apolar) and high (polar) dipole moment 
species ( Sandford et al. . 19881 ). The sublimation temperature of the polar compo¬ 
nent is ~90 K, much higher than that of the apolar component (~20 K), and thus 
the ratio of their column densities is o ften used to trace the ice temperature (Fig. 01 


Tielens et al. 


1991 


Chiaretal 


^ ). 


For the polar CO component, entrapment 
or migration of CO into the abundantly present H 2 O ice seems an evident explana¬ 
tion, but the expected band at 4.647 /im (2152 cm“^) due to dangling OH bonds 
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4.64 4.66 4.68 4.70 


Wavelength [/U.m] 


Figure 4: The CO ice band profile vari es strongly with li ne of sight temperature. 
The Taurus background star Elias 16 ( Chiar et ah . 1995ll and the LYSO EC 90 


( PontoDPidan et ah . 20031 1 trace dust b elow the pure CO sub limation temperature 


of ~ 20 K, while the LYSO L 1489 IRS (Booeert et ah . 2002b 1 and the intermediate 


mass YSO TPSC 78 (jPontoppidan et al.l . l2003l ) have temperatures primarily above 


that. Two laboratory spectra are used to ht the profiles: pure CO corrected for CDE 
grain shapes (green) and CO in polar ices (red). The polar ice is H 20 :C 0 = 4:1 in 
this case, but mixtures with CH 3 OH are likely more appropriate 1 114.21) . The black- 
colored feature on the short-wavelength side is a Gaussian representing a CO ice that 
probably contains CO 2 , for lack of appropriate laboratory spectra. Ice sublimation 
is also evident by the increasing depth of the narrow ro-vibrational lines from gas 
phase CO. 


in interstellar spectra is missing ( Sandford et ah . 1988 I l Mixtures with CH. 3 OH 


rich i ces are more likely, both spectroscopically and chemically 1 118.11 ICuppen et al 


201 il l. High spectral resolution observations have resolved the apolar component 


i nto a blue-shi f ted sho ulder at 4.665 pm and a very narrow feature at 4.673 pm 


( Booeert et ah . 2002blL The latter traces nearly (>90%) pure CO. The detection 


of the solid ^^CO band at 4.78 ^m was pivotal in this assessment ( Boogert et al. 


2002al L The absorption profile of this trace species is insensitive to particle shape 
effects, and very sensitive to the ice composition. For the ^^CO band, grains shaped 
according to the Continuous Distribution of Ellipsoids (CDE) in the Rayleigh limit 
provide the best fit. The composition of the 4.665 pm component is less clear. It 
most likely contains CO 2 , supported by a good correlation of it s strength with that 


of th e CO:C02 component in the 15.2 pm CO 2 bending mode (|PontoDDidan et al. 


20081 1 , but a comparison with l aboratory spectra sho ws that additional species are re¬ 


quired in the mixture as well (jBoogert et al.l . l2002bl l . An alternative explanation by 
absorption of linearly polarized light in a crystalline CO ice phase (the longitudinal 
optical mode; Pontoopidan et al. 2003l l appears less likely. 


Apolar ices: ices 
dominated by 

molecules with low 
dipole moments 
(e.g., CO, N 2 ), 
which have relatively 
low sublimation 
temperatures. 
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Particle shape 
effects: band pro¬ 
file changes caused 
by resonances be¬ 
tween the molecular 
dipoles and the 
particle charge dis¬ 
tribution induced by 
the radiation field 
(iBohren fc Huffman 


Ice distillation: 

sublimation of the 
most volatile species 
in mixtures (e.g., 
CO in CO 2 ice), 
leaving the species 
with higher sublima¬ 
tion temperatures 
behind. 


4.3 The 15.2 /rm CO 2 Band: Segregation and Distillation 

The 15.2 /xm band of CO 2 has proven to be a particularly powerful tracer of the 
composition and processing hist ory of the ices (F i g. El). Building on initial work on 


ISO observations of MYSOs by Gerakines et al. (jl999l i. Pontoppidan et al. ( 2008l i 


use the much larger sample observed by Spitzer toward LYSOs for an empirical 
decomposition of the complex band profiles. All observed bands can be fitted with a 
combination of 5 laboratory spectra, that are fixed in position and shape (Tabled]), 
but whose relative abundances vary as a result of the conditions along the line of 
sight (Fig. Ej): CO 2 (“pure”), C 02 /C 0 ~l (“apolar”), C 02 /C 0 « 1 (“dilute”), 
C 02 /H 20 < 1 (“polar”), and C 02 /CH 30 H< 1 (“shoulder”). The dominance of 
the polar component is strongly supported by the correlation of the CO 2 and H 2 O 
column densities ( ^6.3)) . A subset of YSOs shows distinct peaks at 15.11 and 15.26 
^m, which are attributed to those of pure CO 2 (the symmetry of this linear molecule 
is broken in the pure ice matrix). These are the YSOs with the warmer envelopes, 
and it is thus concluded that pure CO 2 ice is formed by the sublimation of CO from 
apolar ices (“distillation”; 20 K) and at higher temperatures by “segregation” from 
polar ices as t he H 2 O crystallizes and the rearrangement of bonds leads to clusters 
of pure CO 2 (jEhrenfreund et al.lll999l : 30-77 K). A few sources show a third peak 
at 15.35 /xm, which is attributed to bonding of the carbon atoms of CH 3 OH with 


Evolutionary Scenario 
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Figure 5: CO 2 bending mode profiles observed in a range of environments. The col- 
ored areas show the components ( il4.3ll : polar (red), apolar (black), dilute (green), 
pure (yellow), and shoulder (blue). The left panel highlights an evolutionary sce¬ 
nario, showing, from top to bottom, the effects of increasing density and subsequent 
heating. The right panel shows sources that have experienced an extreme heating 
event (HOPS 68 ), have a high CH 3 OH abundance (W 33A), or high ( SST GC 797384) 
and low (LMC 0522 2 2.46) metallicities. The data we r e taken fro m 


(1999), Bergin et al. ( 20^1. Pontoppidan et al.l ( 2008l i. An et al. ( 2011 b 
( 2011 b and Poteet et al. (2013) 
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the oxygen atoms of C 02 _(thn 


of thermal processing (jPartois et al. 


•‘shou Mer” component). This is not an indicator 
I 999 I ). Finally, the “apolar” component is 


probably responsible for the 4.66 pm component of t he CO band 1 ^4.21) as well, and 


indeed their strengths correlate (jPontoppidan et al 


The profiles of the C02-rich components (“apolar” and “pure”) of the strong 
15.2 pm. band are affected by grain shapes in addition to composition and thermal 
history. As for CO, good fits are obtained by CDE shaped grains in the Rayleigh 
limit. T he essential confirmat ion comes from observations of the ^^C02 band at 
4.38 pm ( Boogert et ah . 2000al l. which is insensitive to the grain shapes and can be 
fitted with the same CO 2 components. One striking exception of the CDE model is 
the LYSO HOPS 68, which shows convincing evidence for a population of spherical 
grains, reflecting a n unusual history of sublimation and rapid re-condensation (Pig.EJ 
Poteet et al. 2013 1. 


4.4 The 6.0 and 6.85 pm Features: Complexity Revealed 

The bending mode of H 2 O peaks at 6.0 pm and overlaps with an H 2 O combi¬ 
nation mode that extends to 8 pm. These are much weaker than the 3.0 pm 
stretch mode ltf4.1l). b ecause they benefit less from the hydrogen bonding network 
( van Thiel et al.l . 119571 1. H 2 O absorption being less dominant, the 5-7 pm region has 


proven to be a good tracer of the complex composition of the ices (Pig. ED, despite 
considerable identification issues 1 ^6.11) . 

High spectral resolution observatio ns of MYSOs with IS O revealed variations of the 
peak position of the 6.85 pm feature (jKeane et al.l . l2001bll . which could be attributed 
to the thermal history. Also, a prominent broad absorpt ion overlying bo t h the 
6.0 and 6.85 ^m features was discovered in a few sources (jGibb k. Whittetl . l2002l : 
Gibb et al.l . 120041 1. The 5-7 pm absorption observed toward LYSOs and background 
stars with Spitz erwas decomposed in five components, after the subtraction of pure 
H 2 O ice (Eig. El Boogert et al. 2008l l. essentially confirming the earlier results toward 
MYSOs: Cl (5.84 pm) and C2 (6.18 pm), C3 (6.755 ^m) and C4 (6.943 pm) and 
component C5 peaking at 5.9 pm but stretching over the full 5.8-8 pm range. Each 
component probably represents more than one carrier. The Cl and C3 components 
vary little relative to the H 2 O column density. Component C4 and perhaps C2 
become more prominent at low H 2 O ice abundances, i.e., at temperatures above 
H 2 O sublimation (> 90 K). These features can thus be treated as tracers of extreme 
conditions. Such cases have proven to be rare for LYSOs. Enhancements of the 
C5 component do not appear to correlate with those of C4, and may thus trace 
different processing conditions (e.g., at lower temperatures or by energetic photons or 
particles). Measurements of the strength and shape of this component are inaccurate, 
however, because it is so wide and overlaps with the H 2 O bending mode (Eig. El)- 


4.5 The Lattice Modes: Ice Structure and Emission 


Altho ugh the lattice modes were detected early on (in Orion-KL; lErickson et al 


198 il l, relatively few observations were published due to the limited availability 
of telescope instrumentation at wavelengths above 25 pm. These observations do 
show unequivocally the sensitivity of the modes to the thermal history of the ice 


(Malfait et al. 


1999; 


Demvk et al.l . 120001 : IChiang et al 


Dartois et al 


19981: 


Omont et al.l. Il990l : ISvlvester et al. 


1999 


2 OOIII . The H 2 O feature at 62 pm is particularly 


Ice segregation: 

rearrangement of 
molecular bonds due 
to heating, leading to 
separation of species 
in ice mixtures, 
usually accompanied 
by crystallization. 


ISO: Infrared Space 
Observatory (1995- 
1998), operating 
spectrometers at 
wavelengths of 

2.3-200 pm and re¬ 
solving powers of up 
to R = A/AA ~2000. 


Spitzer: Spitzer 

Space Telescope, 
operating the IRS 
spectrometer (2003- 
2009) at R ~60-128 
(5-14 pm) and 
~ 100 - 600 (10-37 
^m). 
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Figure 6 : Left Panels: Laboratory spectra of species with modes in the 5-8 gm 
region. The sum o f these is compared to the observed 5-8 gm features toward the 
LYSO HH 46 IRS ( Boogert et ah . 2008l b The pure H 2 O contribution is given in red 


and the sum of all laboratory spectra including and excluding the uncertain HCOOH 
identification in green and yellow, respectively. A significant amount of absorption 
is unaccounted for. Right Panels: Ex ample of the decom position of the H 2 O- 
subtracted spectra of the MYSO W 33A ( Keane et ah . 2001b i. the LYSO IC 1396a 
( Reach et al.l . 20091 ). and the LYSO HH 46 in five empirically determined components 
C1-C5, used to investigate source-to-source variations of the band profiles. 


strong in crystalline ices. Laboratory experiments have shown that the bands are 
very se nsitive to the ice composi t ion as well, be c ause impurities affect the long range 
modes (jMoore &: Hudsonl. Il994l : IIoddoIo et al.l . l2014l L In fact, along with the tor¬ 
sional modes at similar wavelengths, they identify complex species more uniquely 
than do the mid-infrared intra-molecular vibrations. Finally, as the lattice modes 
may appear in emission, they show great promise as tools to map the ice properties. 


5 ICES ACROSS ENVIRONMENTS 


Ices are present in a wide range of environments outside of the solar system: en¬ 
velopes and disks of young and evolved stars and Galactic and extragalactic molecu¬ 
lar clouds and cores (Fig. [3l). These environments are characterized by high densities 
(> 10^ cm“^), low temp eratures (< 90 K) and w eak UV radiation fields (Go < 0.07 


at Ay> 1.6 mag; ^8.31 Hollenbach et al. 20091) . The same species are generally 
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detected in all environments, with the exception of the envelopes of evolved stars 
which only show H 2 O ice. The ice abundances vary greatly, however, both across 
and within environment classes. This reflects differences in physical conditions, age, 
and elemental abundances. Observational studies of the different environments are 
briefly reviewed here, and abundances are summarized in ^ 

5.1 Quiescent Dense Clouds and Cores 


The first detections of the 3.0 pm H 2 O ice band toward background stars showed 
that a dense cloud or core environment and a threshold extinction are needed for 


the formation of detectable quantities of H 2 O ices (jWhittet et all . Il983li. In Taurus , 


the ic e band is only present at extinctions above Ay= 3.2 ± 0.1 mag (jWhittet et al. 


200 ll ). which corresponds to an ice formation threshold of Ay= 1.6 mag consider¬ 
ing that background stars tr ace both sides of th e cloud. Similar values are seen 


(Fig. [71) in the Lupu s cloud (IBooeert et ahl. l2013li as well as the dense cores L183 


( Whittet et ah . 2 OI 3 I 1 and IC5146 ( Chiar et al. . 2011). Larger thresholds claimed 


for some environments are not as well established. Toward the Oph cloud, generally 
considered a high radiation environment hindering ice growt h, few background s tars 


were observed and threshold determinations include YSOs (Tanaka et ah. 19901 


Solid CO 2 , which was found to be ubiquitous in quiescent dens e clouds (IWhittet et ahla monolayer has 


1998IL has eff e ctivel y the same formation threshold as H 2 O (iBergin et al.1 . l2005l : formed. 


Whittet et al. 


Ice formation 
threshold: cloud 

extinction depth 
(j4v) indicating the 
onset of rapid ice 
mantle growth, after 


2nn^L showing a close chemical connection between these species. 
Still, their relative abundances vary by up to a factor of three f ll6.3p . Freeze-out 
of the more vol atile CO species o c curs deeper into t he cloud, at ~ 3 mag (Taurus, 


Serpens; Fig. [71 Chiar et al. 1995 I l Kerr et al. ( 1993l f determine a three times larger 
threshold in Ophiuc hus, but along i ndivid ual lines of sight, values are lower, probably 


due to dumpiness (IShuning et al 


depths (9 ± 3 mag; iBoogert et al 


20001) . CH.-^OH ice form s at even greater clou d 

2011blL 


2011; 

Chiar et al. 

2011; 

Whittet et al. 


consistent with the formation requirement of “catastrophic” CO freeze out Bl8.ip . 


rus ( 

Chiar et al. 

. 19961 

Chiar et ah. 

2011). 


to 12% in some isolated dense cores (jBoogert et al 


2011 


A general characteristic is the pri stine nature o f the ices in these environments: 


the dominance of amorphous H 9 O ([Smith et al.l. ll993l i. the large contribution of 


apolar ices to the CO band (Whittet et ah . 19851) . and the absence of segregated 


(heated) CO 2 ices ( Whittet et ah . 20091) . The peak position of the 6.85 pm featur e 


is also consistent with unprocessed ices (iKnez et al 


Knez et ah. 

0 

0 

Boogert et ah. 

2011) 


The absence of the 4.62 pm “XCN” feature ( Whittet et ah . 1985 : Noble et ah . 20131 ) 


is sometimes interpreted as a lack of processing, but, more likely, it is related to 
insufficient CO freeze-out f 116.II and 

5.2 Envelopes of Massive YSOs in the Galactic Disk 


Monolayer: one 

layer of molecules 
covering dust grain 
surfaces, for H 2 O 
corresponding to 
IO 15 niolecules/cm^ 
and an abundance 
Ah ~fewxl0“®. 


A sample of ~20 MYSOs (luminosities above 10^ Lq, masses above 10 Mq) in the 
Galactic Disk has been stud ied extensively in t he ice features, using mos t ly data from 
ISO. W'e refer to reviews bv ICibb et ah! ([20041 ). iBoogert fc Ehrenfreundl (l2004bl ) . and 
Dartois[ (l2005[ ) for historical work. Few new observations of these objects have been 
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Figure 7: Relation between H 2 O (black), CO (green), CO 2 (red), and CH 3 OH (blue) 
ice column densities and extinction for quiescent lines of sight in nearby clouds 
and cores. Upper limits are indicated by open symbols. The lines represent linear 
fits to the detections. For H 2 O, data points are shown for the Taurus (circles and 
solid line) and Lupus IV (asterisks and dash-dot line) clouds, as well as the L183 
(triangles and dashed line) and IC 5146 (squares and dotted line) cores. For CO 
and CO 2 , the data points are for Taurus only, and for CH 3 OH they are for a variet y 
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published since the launch of Spitzer Space Telescope in 2003, but the earlier data 
have been re-analyzed with new methods or in c luded as comparison objects with 
LYSOs in later work ( Pontoppidan et ah . 2008 : Boogert et ah . 2008IL MYSOs in 


the GC region are discussed in il5.5[ 

The evidence for thermal processing of the MYSO ices is ubiquitous in all the 
tracers discussed in IjH This is mostly based on observations along a single line of 
sight toward the central object, but ice maps have confirmed it. Using extended 
scattered light as a background source, the 3.0 gm band is deep toward a disk- 
like structure around the MYSO AFGL 2136, but the ice abundance is suppressed 
relative to the surrounding cloud ( Holbrook &: Temi . IQOsI L probably as a result of 
subli mation. Similar conc lusions were drawn for ice maps toward the MYSO S140 
IRSl ( Harker et ah . IQQtI L 


5.3 Envelopes of Low and Intermediate Mass YSOs 


As the link between quiescent clouds and protoplanetary disks that give birth 
to potentially life-bearing planetary systems, ices in the envelopes of LYSOs are of 
paramount importance. The first (3.0 //m) ice band toward an LYSO was detected 
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by Cohen ( ]j75l ^ toward HL Tau and it was thought to or iginate in an edge-o n disk 
(Cohen, 198,'il ). but a location in the envelope is more likely ( Fiirlan et 2008). This 


demon strates the difficulty of locating ices along the line of sight ([Boogert et al. 
2002cl b Here, surveys of ices most likely located in envelopes around LYSOs are 
reviewed, while disk studies are reviewed in 


Comprehensive surveys of ices in the envelopes of embedded LYSOs were enabled 
by the Spitzer mission (5-20 pm) and sensitive ground-based 3-5 pm instruments. 
The combination of the observations proved critical in the separation of overlap¬ 
ping bands (e.g., in the 5-8 ^m region; 116.1I1 . in obtaining stronger constraints on 
the ice mantle composition (e.g., unraveling the 15.2 pm CO 2 band profile; 114.31) . 
and in obtaining a complete inventory. Samples of ~ 45 objects in nearby clouds 
(Perseus, Serpens, Corona Austral is, Taurus, Oph) an d several cores were studied 


CO9 (PontoDPidan et al.l . l2003l . 120081), CH4 (Qberg et al.l. I 2 OO 8 IL C HqQH and NH3 


in th e 5-8 pm and H^O features (iBoogert et al.l. 120081) and th e bands of CO and 


( Bottinelli et ah . 201C ). and “XCN” ( van Broekhuizen et ah . 20051 ) . A survey of all 


i ce absorptions i n the 5-20 pm range focused on LYSOs in Taurus was done as well 


( Zasowski et ah . 20091 ) 


These studies consider mostly Class 0 and I LYSOs. As for MYSOs, variations in 
abundances and processing history were observed. Thermal processing was found to 
play an important role, despite the lower stellar luminosities. Aided by episodic heat¬ 
ing events, sublimation of the most volatile C O and CO 2 comp onents was observed 
for even the lowest luminosity YSOs (<0.1 L^T^ uKim et al.ll 2012 l ). Large variations of 
the OCN“ abundance (4.64 um ), even for sources close toget her on the sky, suggest 
special conditions in envelopes ( van Broekhuizen et ah . 20051 ). They may be related 
to past CO freeze-out and possibly heating events (^. 

LYSOs in clustered and high radiation environments of nearby ma ssive stars ma y 
be more representative of the early solar system than isolated YSOs ( Adams . 201ol ). 
An example is the IC 1396 globule irradiated b y a nearby 0-star. It shows 6.85 pm 
features peaking at relatively long wavelengths ( Reach et ah . 20091 ) . a high temper¬ 
ature effect that is rarely observed toward isolated LYSOs (' 114.41) . 

Finally, observations of the envelopes of intermediate mass YSOs (masses of 2 - 
10 Mq and luminosities of a few to 10^ Lq) would enable further studies of the 
dependency of ice characteristics on stellar mass, but the sample sizes are small. In 
a Spitzer survey of 14 embedded intermediate mass YSOs distributed on the sky, 
Pitann et al.l (| 201 ll ) showed that the common ice features are present and that the 
15.2 pm CO 2 band shows no signs of segregation. Thermal processing is, however, 
evident by r educed H^O and CO ice abundances in intermediate mass YSOs in the 

One of these sources has high CH 3 OH and OCN“ 


Vela cloud (jThi et al 


> 00 ^ ). 


abundances, but this is not necessarily a sign of thermal processing f ^ 6 . 1 l and l 8 .ip . 


5.4 Disks Around Low and Intermediate Mass YSOs 


Direct observations of ices in circumstellar disks trace the volatile building blocks 
of planet and comet formation. These observations are challenging because the rel¬ 
atively flat disk geometry requires a specific orientation for the 2-20 pm vibrational 
modes to be seen in absorption. Great improvements in the sensitivity and spa¬ 
tial resolution of infrared observations, as well as in disk models, has led to much 
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progress, however. A recurring conclusion is that ices in disks are strongly affected 
by thermally processing. 

For flared disks, the observed ice features are most pro minent if the disk incl i natio n 


is comparable to the opening angle of about 70 degrees (jPontoDPidan et al.l. 120051 1 


Lower angles trace less dense outer disk regions and the observed ice features are also 
diluted as the size of background scattered light region exceeds that of the absorbing 
region. At higher angles, the high disk mid-plane densities increase the dust optical 
depths. As a result, ice observations are biased toward processed ices in the disk 
surface l ayers. This is indeed obser ved in the edge-on disk source CRBR 2422.8-3423 


in Oph (jPontoppidan et al. 


200 


^), as it is one of the few low mass sources with a 


shifted 6.85 gm feature indicative of high temperature processing (<14.41]. 


R ather than us i ng the light from the central object, iTerada &: Tokunaeal (j2012al l 
and Terada et al. ( 2012bl l use the extended background of the Orion nebula to trace 
the ices. This work confirms the optimal disk inclination angle of 65-75 degrees 
needed to detect ice features. One source shows an unprecedented H 2 O ice band 
profile (Fig. ED with a strong contribution from large (~ 0.8 gm) grains rich in crys¬ 
talline ices that may have surfaced from the lower disk regions where coagulation is 
more efficient. The same authors have resolved a (different) edge-on disk against the 
extended background and find spatial variation s in band depth but not shape. The 
depth does not vary on a time scale of 4 years ( Terada et al. . 2012bl l. implying that 
the crystalline ices are extended over at least 3 AU at distances b eyond 30 AU. This 
distan ce estimate comes from spatially resolved spectroscopy by fschegerer &: Wolf 
( 2 OI 0 I 1 toward the T Tau object YLW 16A. 

High spatial resolution observations have indeed been very powerful in locating 
ices in disks. Coronographic adaptive optics (AO) images in narrow band filters of 
the Herbig Ae star HD 142527 reveal a spatially exten ded 3 /rm absorption attributed 
to scattering by icy ~ 1 gm sized grains in the disk ( Honda et ah . 2009l l. Also, AO 
spectra separated the HK Tau and HV Tau binary and tertiary systems, revealing 
deep (amorphous) ice bands in the HK Tau B and HV Tau C edge-on disks. The 
depth toward HV Tau C varies on a time scale of 2.3 years, which co uld be caused 


by a 1.4 AU sized structure at a distance of 100 AU from the star (jTerada et al. 

2OO7I L 


A particularly successful technique to trace ices in circumstellar disks, that is 
less dependent on the disk orientation, is spectroscopy of the far-infrared ice lattice 
modes ( l)4.5p . In passively heated, flared disks, the lattice modes of H 2 O ice at 44 
and 62 u m are excited in su per-heated optically thin disk surfaces at radii beyond 
100 AU (jChiang et al.l . l200lh . In this layer, the ice condensation radius (not to be 
confused with the “snow line” in the disk mid-plane) scales with the 3rd power of 
the stellar temperature, which explains the detection of the lattice modes in ISO 
spectra of the coolest members of their sample of Herbig AeBe and T Tauri stars. 

The ISO spectrometers detected more of the lattice modes, mostly from crystalline 
H 2 O, in bright disks around i ntermediate mass stars: the isolated H erbig Ae stars HD 


100546 ( Malfait et al. 


OOsI ) and HD 1 42527 (Malfait et al. 


Ae star HD 163296, but not AB Aur (|van den Ancker et al. 


I 999 IL and the Herbig 
2 OO 0 I 1 . Observations 


with the Herschel/PACS spectrometer are limited by calibration uncertainties and 
incomplete coverage of the strongest (44 //m) lattice mode. Detections claimed for 
a few T Tauri stars indicate a high d egree of crystall i nity c reated by past heating 


events, such as planetesimal collisions (jMcClure et al.l . l2015l l. 
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Finally, sensitive, but low spatial resolution Akari observations of a sample of edge- 
on disk sources confirm the optimal inclination that is needed t o detect ice features , 
but the spectral resolution is too low to study processing effects (jAikawa et alJ . l 2012 ll . 
Their 2-5 pm spectra reveal solid H 2 O, CO, OCN“, OCS, and possibly HDO. 

5.5 Galactic Center Region 

The conditions in the GC region differ significantly from those elsewhere in the 
Galaxy (elemental abundances, radiation fields, dynamic time scales) and thus offer a 
unique perspective on ice formation and evolution processes. Studies of ice features in 
the GC region have primarily targeted the cluster of infrared sources (GC IRS) within 
a parsec from SgrA*, and MYSOs (candidates) throughout an extended (~ 50x200 
pc) region in the Central Molecular Zone (CMZ). 

Most of the infrared sources in the central parsec of the Galaxy s how ice absorption 


Akari: a Japanese 
space telescope, op¬ 
erating the InfraRed 
Camera (2006-2011) 
equipped with, 

among others, a 
prism at R ~20-100 
(2.5-5 pm). 


features. The depths of the 3.0 um H 2 O and 4.67 pm CO bands (jMcFadzean et al 


198S : 


2001 


Chiar et al 


2002 : Moultaka et ah . 2005) a nd the 4.62 pm XCN feature ( Moneti et al 


Chiar et ah . I2OO2I : iMoultaka et al.l . bood ) vary with factors of ~4 between 


these sources. As much as 60-80% of the H 2 O ices may be local to the inner pc 
of the Galaxy, pos s ibly in dense cl u mps in the circum-nuclear ring or mini-spiral 


(jChiar et al. 


2002 


certain for CO ices (| Moult aka et al 


Moultaka et al.l . l200.^ 1. The foreground contamination is less 
20091 b but in any case the total column relative 


to H 2 O is small (< 10%). While not much CO has frozen out, there is little evidence 


for thermal processing . The CO profile has a strong apolar component (jMoneti et al. 


2001 


Moultaka et al.l . l2005l ) . and the 15.2 pm CO 2 bending mode, observed in the 


ISO beam of ~ 20 ” encompassing a numb er of IRS sources, shows a lack of thermal 


processing as well (iGerakines et al 


1999lb Other detected ice sp ecies include NH^ 


and CH 4 , while solid CH 3 OH is notably under-abundant (< 5%; IChiar et al.ll 200 fll . 
2002l b The detection of the stretch mode of NH 3 at 2.96 pm in ISO and high spa¬ 
tial resolution ground-based observations, supported by the presence of NH 3 .H 2 O 
hydrates at 3.2-3.7 pm and by the 6.2 pm bend mode, is unusual and implies a large 
NH 3 abundance of up to 30% relative to solid H 2 O. An N-rich enviro nment seems 
also reflected in a high OCN~ abundance of 3% toward GC IRS 19 (IChiar et al 


2002 , recalculated using the band strength in 


Broekhuizen et al 


200 , 5l b 


On a larger scale, ices were detected in numerous sources throughout the CMZ. 
In a sample of 107 targets selected using 2MASS and Spitzer 1-8 pm broad band 
photometry, 3 5 were identified as MYSOs (or candidates; 8-23 Mq) following Spitzer 
spectroscopy (|An et al.l. l201llb They are located throughout the CMZ, although 
there is clustering near the Sgr B molecular complex. All objects show ice features 
in the 5-7 pm region, as well as the 15.2 pm CO 2 bending mode. The latter was used 
as a criterion in the classification of the sources as MYSOs: the presence of a long 
wavelength wing tracing the CH 30 H:C 02 complex was required (Fig. [5]). It is based 
on the fact that until then this wing had only been detected towards Galactic Disk 
MYSOs with high GH 3 OH abundances 1 H4.3D . Further work is needed to confirm 
the nature of these sources and the location of the ices in the MYSO envelopes or 
in the dense cloud medium in the GMZ. Also, spectra at wavelengths shorter than 5 
^m have not yet been obtained, thus direct measurements of the H 2 O, GH 3 OH, and 
CO abundances toward the CMZ MYSOs are lacking. 


























































































22 


Boogert, Gerakines, Sz Whittet 


5.6 Shells Surrounding Evolved Stars 

Certain classes of evolved stars show H 2 O ice bands and offer a view on ice proper¬ 
ties from a completely different perspective compared to that of the ISM and YSOs. 
In general, these are carbon-poor AGB stars with very high mass loss rates, showing 
OH maser emission (OH/IR stars), or the slightly more evolved post-AGB stars that 
still have dense, but expanded and cooled envelopes. In short, the expanding dense 
shells of gas expelled by evolved stars cool rapidly. At some point, H 2 O is formed in 
the warm gas, and as the shell cools it condenses onto the refractory grains formed in 
an earlier, hot ter phase. The 3.0 um band is narrower in OH/IR stars compared to 


YSOs (Fig. [3l Gillett fc Soifei 197d l and is th ought to reflect the purity of the H 2 O 


ice, its crystallinity, and the small grain size (ISmith et al.l. Il988lb The 6.0 bending 


mode ( Soifer et ah . 198lh . the 11 gm libration mode ( Sylvester et ah . 1999l i and the 
lattice modes at 44 and 62 gm were detected as well. The lattice mode was first 
detected by anomalously str ong IRAS 60 /rm broad band emission, half of which 
can be attributed to H 2 O ice ( Forveille et ah . 1987 : Omont et ah . 199ol h in the post- 
AGB star Frosty Leo. The lattice modes strongly confirm the crystallinity of the ice 


(jOmont et al.l . 11990 : ISvlvester et al.l.ll999l:lDemvk et al.1. 120001). No other ice spe cies 
were detected toward evolved stars ( Gillett &: Soifei . 1976 : Sylvester et ah . 1999l i. in 


part due to the oxygen rich environment in which they were formed. 


5.7 Massive YSOs in the Magellanic Clonds 


Observations of ices in the Magellanic Glouds offer unique tests of ice and molecule 
formation and evolution models, due to the lower metallicities and correspondingly 
lower dust content and higher UV fields. Overall, differences in the relative ice 
abundances and higher levels of processing have been found relative to Galactic ices. 

The GO 2 /H 2 O column density ratios are a factor of two larger in the envelopes of 


MYS Os (10^ — 10® L m) in the Large Magell a nic Gloud (LMG: 33T1%: Seale et al 


2011 


using data from IShimonishi et al.l l2008l: lOliveira et al.ll2009l : IShimonishi et al 


201 oi l compared to the Galaxy (17±3%; 


Gerakines et al 


15.2 gm GO 2 band is very well studied 1 ^6.31 


Seale et al 


lOodh The shape of th e 

2009l i. 


2011 


Oliveira et al 


In the majority of sightlines, GO 2 ices have experienced th ermal processing. GO 
ice is detected in ground-based s pectr a (|01iveira et al.l . 120111 1 and in low resolution 
Akari spectra ( Shimonishi et ah . 201ol l. indicating a low abundance (< 15% relative 
to H 2 O). The 5-7 gm ice features, including a clear 6.85 gm feature in at least 
one object, were detected ( Oliveira et al.l. 2009) although t heir analysis is severely 


hampered by overlapping PAH emission ( Seale et ah . 2 OO 9 I 1 . The XGN feature at 


4.62 ^m m ay have been detected in a few Akari sources but suffers from low spectral 


resolution ( Shimonishi et ah . 201ol l. This is the case for the 3.53 gm. band of solid 


GH 3 OH as well. 

H 2 O and GO 2 ices were also detected toward 14 MYSOs in the S mall Magellanic 


Glou d (SMG) in a combined Spitzer and ground-based 3-4 gm survey (jOliveira et al. 


2013l i. The relative column densities are compared to other environments in 
Heating effects play a role here as well, judgi ng the presence of the 62 gm lattice 
mode of crystalline H 2 O ice in Spitzer spectra (jvan Loon et al.l . 120101 ). The 6.0 and 

6.8 ^m features were detected, but GO ice is absent at levels below those in the LMG 
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(although only four lines of sight were observed; Oliveira et al. 201 ll ) 


5.8 Other Galaxies 

Ice features are commonly seen in the infrared spectra of external galaxies with 
centers dominated by starbursts or dust-enshrouded Active Galact ic Nuclei (AGN). 
They are often used to determine the nature of the emission sources. Ilmanishi &: Malonev 
( 2003 ! ) show that peak optical depths of the 3.0 pm band larger than ~ 0.3 must 
be caused by compact AGN emission surrounded by dust clouds containing the ices. 
Weaker ice features are due to starburst regions in which many MYSOs are mixed 
with clouds containing ice coated dust grains within the same beam. Local ULIRGS 
and most LINER-type galaxies show deep ice bands, and their centers are thus dust- 
enshrouded AGN. In contr ast, Seyfert-type gala xies do not show ice bands at all, 


lacking such dusty centers ( Imanishi et al. . 20od L 


Besides H 2 O ic e (3.0 and 6.0 um), solid C O 2 , “XCN”, GO and CH 4 were detected 


in these galaxies (jSDOon et al.l . I 2 OOOI . l2004l ). The spectra often show a combination 
of PAH emis sion, and absorptio ns from diffuse dust and ices, much like the center of 
the Galaxy ( Sturm et ah . 200 oI l complicating the analysis of the ice features. 

The ice distribution was mapped in several nearby galaxies, and enhanced relative 
abundances are often attributed to production in energetic radiation fields. Toward 
a ring encircling the starburst nucleus of NGC 4945, the 4.62 pm XGN feature 
is stronger than that of CO, similar to t hat in the G a lactic MYSO W33A. The 
depths do not vary much within the ring. Spoon et al. ( 2 OO 3 I 1 attribute it to UV- 
processed ices (X-ray fields are likely too weak). The Akari satellite mapped the 
ice distributions in other starburst-dominated galaxies. In M82, variations of the 
CO 2 /H 2 O abundance ratio correlate with the U V field traced by Bra emission, 
possibly indicating the energetic production of CO 2 ( Yamagishi et ah . 20I3I ). In view 
of the difficulty of identifying energetically-produced ices in Galactic environments 
these claims for nearby galaxies must be regarded as te ntative. _^ 

2nn9lL 


The current redshift record for ice feature detections is 2.28 (jSaiina et al. 


At such high redshift, both the 3 and 6 pm ice bands are within the 5piteer spectrom¬ 
eter wavelength range. Overall, no differences in the ice composition or abundance 
were found as a function of redshift, but the data are limited. 

5.9 Solar System Objects 


There are many ice-rich environments in the Solar System, including the satellites 
of the gas- and ice-giant planets, trans-Neptunian Objects (TNOs), and comets. 
For this review of interstellar and circumstellar ices, cometary ices are of particular 
interest because of their possible link with ices in the ISM. Ice abundances are the 
main observables to test this link. While present day comets originate from distinct 
heliocentric distances, formation regions are less clear as a result of radial mixing. For 
example, organic-enriched comets are present in both Kuiper Belt and Oort Gloud 
comets. A certain population of comets may even be captured from sibling star 
systems formed from the same cluster as th e Sun. Gometary ice abund ances will not 
be reviewed here, but rather the review bv iMumma &: GharnlevI (|20Ilh will be used 
as a reference. Remarkable differences and similarities with respect to interstellar 
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MRN: number den¬ 
sity distribution of 
grain radii derived 
from diffuse cloud 
observ ations, n(a) oc 
a~^'^ ( Mathis et al. . 

1973). 

Upper and lower 
quartile abun¬ 
dances: median of 
the abundances re¬ 
spectively larger and 
less than the median 
of all abundances. 


ices are observed ( (16.51) . 

6 COLUMN DENSITIES AND ABUNDANCES 


Ice column densities are generally determined by dividing the integrated optical 
depth by the integrated band strength A measured in laboratory experiments (^: 


A = y Tudv/A [cm 


( 1 ) 


For saturated bands, and in particular for unresolved ones such as those observed 
with the low resolution mode of the Akari satellite {R ^ 20), a curve of growt h 


analysis must be performed to derive the column density (jShimonishi et al 


For LYSOs with disks, the size of the background region of scattered light, and thus 
the observed ice column in the foreground, may be a strong function of wavelength 
( PontoDPidan et ah . 2005lh Different modes of the same species would then yield 


different column densities, and in fact their ratio would trace the disk inclination. 
So far this has not been observationally demonstrated, however. Similar scenarios 
were suggested to explain discrepan cies between the 3 a nd 6 gm absorption band 
depths towards MYSOs due to H 2 O ( Keane et ah . 2001bl i. but this also still remains 
to be proven. For example, toward the MYSO S140 IRSl, the column densities 
derived for al l CO? b ands over the 2-15 gm wavelength range are in good agreement 
([Keane et al.l . l2001al h Finally, for external galaxies such as ULIRGs and starburst 
galaxies, a single spatial resolution element may contain multiple continuum sources 
(e.g., YSOs, an AGN), and depending on the geometry of those and the dense clouds 
conta ining ices, the ice band depths may be quite different ( (14.11 llmanishi fc Malonev 

20o3). 

Equation [1] only applies to pure absorption. For features caused by scattering 
(^, the ice column traced depends strongly on the grain size, and radiative transfer 
involv ing grain size distribut i ons and appropriate optical co nstants must be calcu¬ 
lated ([Pendleton et al.l . [l990|; iDartois &: d’Hendecourtl . 120011 1 . Similarly, the lattice 
modes (> 25 gm) may appear in emission and radiative transfer models are re- 
quired to evaluate co lumn densities such as demonstrated f or disks surrounding YSOs 


1996 


( Chiang et al.l . 200ll l and envelopes surrounding MYSOs ( Dartois et ah . 

H 2 O ice was detected at column densities ranging from 10^^ to 10^® cm The 
lower limit corresponds to a peak optical depth of the 3.0 gm band of 0.06, which is 
rather high because of uncertainties in the continuum determination. It corresponds 
to dust extinctions close to the ice formation threshold ( ^5.ip and thus to a few 
monolayers of ice on each grain in th e MRN grain size distribution with a lower 
limit of 20 A ( Hollenbach et ah . 20091 ) . H 2 O might be present on grains below the 
ice formation threshold, but, besides the low column, the observational signature is 
weakened by the much lower band strength of isolated water compared to that of 
bulk ice ( (I4.ip . 

Ice abu ndances relative to H 2 O, AhjOj were determined following the statistical 
method in Oberg et al.l ( 201ll ). but using somewhat different samples and including 
minor ice species. The values for the Galactic sources are listed in Table [2j For 
samples larger than a few, the median and lower and upper quartiles of the abundance 
distributions are given, as well as the median abundances taking into account upper 
limits. Also listed are the full ranges of the observed abundances. 
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Cometary abundances derived from measured gas phase production rates were 
taken from Mumma fc Charnlev ( 2f)lll ). supplemented by measurements of CO 2 
abundances with th e Akari satellite (within 2.7 AU, the H 2 O sublimation radius; 


Ootsubo et al. 2012I '). Considerable variations are observed within and across envi¬ 


ronments f li6.3p . The insecurely identified species are discussed in ^6.11 and upper 
limits for other astrochemically relevant species in 














Table 2; Ice Abundances 


ts3 

ai 


Species 



[%] 



Y9 [10-5] 



MYSOs 

LYSOs 

BG Stars® 

Comets 

MYSOs 

LYSOs 

BG Stars® 

Securely identified species'^: 






H2O" 

100 

100 

100 

100 

3150 

3855 (42) 

401? (39) 






12-57 

14-80 

(< 9)-62 

CO® 

(7) 

2135 (18) 

25i 


2.611 (1-9) 

9.6^8 (8.1) 

1220 


3-26 

(<3)-85 

9-67 

0.4-30 

(< 0.4)-12.8 

(< 1.2)-26 

3-21 

CO2" 

i9fi 

2811 

2639 

1524 

-‘■OlO 

3.7^ 

11.8i°o 

13.238 


11-27 

12-50 

14-43 

4-30 

1.8-15.6 

2.4-25 

5.2-26 

CH3OH 

Of (5) 

6P (5) 

8^9 (6) 


3.7139 (1.7) 

3.31:9 (2.3) 

5.211 (2.4) 


(< 3)-31 

(< l)-25 

(<1)-12 

0.2-7 

(< 0.4)-16.6 

(< 0.2)-15 

(< 0.6)-6.6 

NH3 


6l (4) 




3.61:1 (2.6) 



~7f 

3-10 

< 7 

0.2-1.4 


(< 0.4)-6.4 

< 4 

CH4 


4.5i (3) 




2.33:1 (1.4) 



1-3 

1-11 

< 3 

0.4-1.6 

0.4-1.8 

(< 0.2)-5.6 

<1.2 

Likely identified species® 







H2CO 

~2-7 

~6 


0.11-1.0 




OCN- 

0 60-7 

0.6[]:| (0.4) 



0.41:1 

0.41:1 (0.06) 



0.1-1.9 

(< 0.1)-1.1 

<0.5 


0.2-0.6 

(<0.04)-0.4 

<0.16 

OCS 

0.03-0.16 

<1.6 

<0.22 

0.1-0.4 

(< 0.02)-0.06 

<3.2 

<0.12 


Possibly identified species’^: 


HCOOff 4| (3) 

(< 0.5)-6 (< 0.5)-4 <2 

CHgCHaOff ~Xh2o(HCOOH) 


0.06-0.14 


2.12| (0.9) 2.4f| (0.9) 

(< 0.08)-2.6 (< 0.4)-7.8 

~ a:h(hcooh) 


<4 
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Species 

MYSOs 

[%] 

LYSOs BG Stars'^ Comets 

MYSOs 

Ab [ 10 - 6 ] 
LYSOs 

BG Stars^ 

HCOO-j 

o.5g:I (0.5) 




O.ISO;?! (0.14) 




0.3-1.0 

~0.4 

< 0.1 


(< 0.10)-0.42 

~ 0.2 

< 0.8 

CH 3 CHOj 

Ah 2 o(HCOO-)x 11 



Ah(HCOO-) 

xll 


NH+ 

11 13 

Ills 

sii 

06 


4 lS-4 

4.6i:| 

q q4.9 


9-34 

4-25 

4-13 


1.4-6.0 

0 . 8-12 

1.9-9.6 

SO 2 

(< 0.9)-1.4 

~ 0.2 


0.2 

<0.4 

~0.08 


PAH^ 

~8 








^ Abundances relative to the H 2 O ice column density. For each molecule, the first row gives the median 
and lower and upper quartile values of the detections, and in brackets the median including upper limits. 
The second row gives the full range of abundances. If only one row is given, it shows the full range. See 
Table [Hand ^for references. 

As in the Xn^o column, but for abundances relative to Ah- 
^ Quiescent clouds and cores. 

^ Identified based on multiple modes or isotopologues in high quality spectra. 

® The isotopologues and ^^C 02 were also securely identified, and HDO was possibly identified 


^ Values derived from the simultaneous analysis of multiple features ( Dartois et al.l . 120021 ) . 

® Identification based on a single absorption feature and the profile matches laboratory spectra. 

^ Identification based on a single absorption feature and no convincing match to laboratory spectra is 
available. 


^ HC OOH and/or CH 3 CH 2 OH may be carriers of the 7.24 //m ice band (jSchutte et al.l . ll999l : l6berg et al. 

201 ih . ^_ 


199 


3 ). 


J HCOO“ and/or CH 3 CHO may be carriers of the 7.41 um ice band (ISchutte et al. 

Using column density of C-C bonds derived in iHardegree-Ullman et al.l ( 20141 ) and assuming each PAH 
species contains 50 carbon atoms. 


ts2 
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Ice abundances relative to elemental hydrogen, Xu, are generally calculated by 
determining A^h = A^(HI) + 2A^(H2) from scaling relations with the silicate fea¬ 
ture (r 9 . 7 ) or the extinction in the near-infrared (^k) or visu al (Ay). These abun - 


dances are also given in Table [H assuming r g 7 = 0.26^k ( Boogert et ah . 2 OI 3 II 


and A^h = 1-54 x 10 ^^Tk ( Vuong et ah . 2003l b Another way to express ice abun¬ 


dances, of particular interest to planetary disk models, is the ice-to-rock mass ratio. 
Defining “rock” as the silicate dust comp onent, the highest observ ed ice abundances 


correspond to ice-to-rock ratios of ~1.5 (jPontoppidan et al.l . l2014l h 


Finally, for extragalactic sources only CO, CO 2 , and H 2 O ices were observed in 

0I8 


4)-19) for the LMC ( 

Shimonishi et ah. 2010: 

Oliveira et ahl. 

2011 

), and < 3 for the 

SMC ( 

Oliveira et ah. 

2011). AnoofCO?) is 17o^ (6-41) for the SMC (Oliveira et ah. 


2013) and 36 oq f 1(3-46) for the LMC K ( Shimonishi et ah . 2nid : Seale et al 


2011 


Oliveira et ah . 20091 : Shimonishi et ah . 20081 ). Abundances relative to hydrogen were 


not published. The Magellanic Cloud abundances are compared with Galactic envi¬ 
ronments in 


6.1 Identification Issues 


The origin of a number o f ice absorption feature s is u ncertain (Table [T]). As 


PAHs: 

Aromatic 

carbon 


opposed to gas phase features ( Herbst &: van Dishoeck . 20091) 


Polycyclic 

Hydro¬ 

molecules, 


the identification of ice 
features is complicated by the dependence of the peak position on the ice environment 
and IH). Therefore, Table [2] is divided into categories of “securely”, “likely”, 
and “possibly” identified carriers. For securely identified species, more than one 
vibrational mode was detected and the depth and profile were successfully fitted with 
models and laboratory spectra. Although CO ice has only one mode (4.67 /un), the 
^^CO isotopologue at 4.78 gm was detected as well. Species in the “likely identified” 
category have only one securely detected absorption band, but an excellent fit to 
the profile and position was demonstrated and the abundance is reasonably well 
understood. Species in the “possibly identified” category lack a convincing profile 
fit or would have an unreasonably high abundances (e.g., compared to gas phase 
measurements). SO 2 is in this category as well, because, although the 7.63 //m 
feature was well fitted by SO 2 in CH 3 OH ices, it is so far detected in only one line 
of sight. In the remainder of this section some identification issues are highlighted. 

The 3.25 gm (PAH?) feature. The 3.25 gm absorption feature depth, mea- 


responsible for the 
ubiquitous infrared 
emission features, 
and, in dense en¬ 
vironments, likely 
also frozen out in ice 
mantles. 



erogeneity of PAH spec ies, including size distributions and ionization within the ice 

An alternative propose d carrier is NH^, which is also less 


(Bouwman et al 


2011 


volatile than H 2 O ices (jSchutte fc Khannal . l2003l ) . and whose spectrum also matches 


the 6.85 /rm feature (see below). 
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The 3.47 /xm (NH 3 ) feature. This feature, not to be confused with the 3.40 pm 



overall long-wavelength wing of the 3.0 pm band and other NH 3 features throughout 
the infrared spectra of YSOs. 


The 4.62 pm (“XCN”) feature. The 4.62 feature was attribut e d to a 


carrier with a C=N (triple bond) group since its first detection (jLacv et al.l . Il984li 


Laboratory experiments in which CO-bearin g ices containing NH 3 were photolyzed 
reproduced this feature ( d’Hendecourt et ah . 19861) and pointed to the C=N bond 
of a small molecule. Hence, the observed absorption became known as “the XCN 
feature”, but m any laboratory stu dies have pointed to the true carrier as the ionic 
species OCN“ (jDemvk et al.lll998l and references therein). A good correlation be¬ 
tween the gas phase HNCO, presumabl y photodesorbed at cloud edges, and solid 


state OCN abundances supports this (jOberg et al.l. l2009al) . Later it was shown 


that the interstellar absorptio n is a composite of features peaking at 4.60 and 4.64 
pm ( PontODPidan et ah . 20031 ). Only the peak position and profile of the 4.64 pm 
component are in agreemen t with the OCN~ species, mo st likely in polar ices that 
have possibly been heated (Ivan Broekhuizen et al.l . l2005l ). It thus appears that its 
strength relative to apolar ices is an indicator of thermal processing, but only for 
environments that have experienced the “catastroph ic” CO freeze out s tage in which 
HNCO and subsequently OCN“ were formed l^S.lMOberg et al.ll201ll) . The nature 
of the 4.60 ^m compon ent is less clea r . iPontoppidan et al.l (j2003l ) suggest CO in a dif¬ 
ferent binding site, and lOberg et al.l (j201ll ) suggest O CN~ in apolar ices alth o ugh n o 
good fit was obtained for any mixture considered by Ivan Broekhuizen et al.l ()2005l ) . 


The 6.0 and 6.85 pm (salts, organic residue, PAH, NH^?) features. 

Securely identified species do not fully explain the observed prominent 6.0 and 6.85 
pm absorption features in most sightlines (Fig. [ 6 ]). Dilution of H 2 O in other species, 
such as the abundant CO 9 . enhances th e strength of the H 2 O bending mode relative 
to the 3.0 pm band ( Knez et ah . 20051 ). but this is generally insufficient to explain 
all 6.0 pm absorption. HCOOH could be an important contributor but its identi¬ 
fication in the ices is uncertain (see below). The carrier of the 6.85 pm feature is 
not related to the bending modes of aliphatic hydrocarbons seen at similar wave¬ 
lengths in the diffuse ISM due to the absence of the corresponding stretching modes 
at 3.40 pm. The CO/ ion in minerals (carbonates, e.g., MgC 03 or CaC 03 ) is also 
excluded due to the absence of features at longer wavelengths in the spectra of YSOs 
(jSchutte et al.l . ll996l) and the excelle nt correlation between the dept hs of the 6.85 pm 
absorption and known ice features (jTielens &: Allamandolal . Il987al ) . The latter ar¬ 


gume nt also excludes other refractory carriers. The NH^ ion in salts ( Knacke et al. 
19821 ). is considered the most likely carrier, but it lacks a co nvincing prof i le fit . 


While the observed strengths correlates very well with H 2 O ice (|Pontoppidanl . l2006l : 
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Boogert et alJ . 120081 ) . the laboratory profile is too broad in H 2 O ices (jGalvez et al.l . 


2010l ). After H 2 O ice sublimation, a better fit is obtained, and at even higher tem¬ 
perature the peak position shifts to longer wavelengths, consistent with observations 
of YSOs (enhanced C4/C3 component ratios; 114.4I] . Finally, a very broad absorp¬ 
tion o bserved toward a number of YSOs (C5 compon ent: 114.41) has been ascr ibed to 
salts ( Schutte &: Khanna . 20031) and organic residue ( Gibb &: Whittet . 2002 ). both 
of which are energetic and/or thermal processing products, but these identifications 
need independent spectroscopic confirmation. 


The 7.24 /rm (HCOOH?) feature. This feature was attributed to the G 
H deformation mode in HGOOH, and specific mixtures cont aining both H 2 O an d 


GH 3 OH are needed to match the peak position and prohle (ISchutte et al. 


199! 


3 ). 


The depth is consis tent with the much stronger G=0 stretch mode at 5.85 gm 
( Schutte et ah . 19961) . but its identification was put into question because of issues 


with the baseline in the laboratory spectra ( Boogert et ah . 2008 : Oberg et ah . 201 ll) 


and an abundance that is a f actor of 10 '^ larger com p ared to the gas pha se in hot 


cores is also not understood (jBisschop et al.l . 120071 ). lOberg et al.l (j201ll ) propose 


pure ethanol instea ,d. Other species te nd to have sli ghtly blue-shifted peak positions 
(HGONH 2 , G 5 H 12 ; Schutte et al. 1999 : NH 2 GH 2 OH; Bossa et ahlbood ) and are listed 
in Table El 


6.2 Upper Limits 


Table [3] shows abundance upper limits of astrophysically relevant species. These 
were determined by comparing laboratory spectra with the spectra of YSOs and 
background stars, taking into account position and profile dependencies on the ice 
composition. For symmetric species (in particular O 2 , N 2 , and H 2 ), interactions 
in the ice matrix induce a dipole moment, and the band strength is a very strong 
function of molecular environment. These make the derived abundance upper limits 
uncertain. Indirect measurem ents of the Oy and N 2 abund ances in mixtures with 
GO by the effect on the ^^GO ( Elsila et ah . 1997 ) and ^^GO ( Boogert et ah . 2002al) 
band profiles lower the upper limits for N 2 by a factor of ~ 5, but do not further 
constrain those of O 2 • The case of H 2 refers to energe tically processed ices that are 
able to retain the produced H 2 (|Sandford et al.l. Il993l ). In pure form H 2 sublimates 
at 3 K, well below the temperature of dense clouds. An upper limit rather than 
the detection claimed in ISandford et al.l (jl993l ) is quoted because of confusion with 
photospheric absorption lines at 2.415 gm. 
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Table 3: Abundance upper limits of astrochemically relevant 
ices and corresponding detections in comets 


Species 

VpaO 

% 

Vh 

10-® 

Environment®' (reference)*’ 

^HaO (comet) 

% 

N 2 

< 0.2 - 60 

< 0.1 - 28 

Taurus cloud (1, 2, 3) 


O 2 

< 39 

< 60 

LYSO R CrA IRS2 (4) 



< 15 

< 30 

MYSO NGC 7538 IRS9 (4) 


H 2 

< 68 

< 14 

LYSO WL 5 (5) 


H 2 S 

<0.3-1 

< 0.04-0.12 

MYSOs (6) 

0.12-1.4 


<1-3 

< 0.6 - 1.6 

Taurus cloud (6) 


H 2 O 2 

<2-17 

<0.6-8 

YSOs, Taurus cloud (7) 


C 2 H 2 

<1-10 

<0.4-4 

MYSOs (8) 

0.1-0.5 

C 2 H 6 

< 0.3 

< 0.14 

MYSO NGC 7539 IRS9 (8) 

0.1-2 

C 5 H 12 

< 15 

< 10 

MYSO W 33A (12) 


C 3 O 2 

< 5 

< 2 

YSOs (9) 


N 2 H 4 , N 2 H+ 

< 10 

< 4 

MYSOs ( 8 ) 


HNCO 

< 0.3-0.7 

< 0.10-0.24 

MYSOs (11) 

0 . 02 - 0.1 

HCONH 2 

< 1.5 

< 1 

MYSO W 33A (12) 

0.002 

NH 2 CH 2 OH 

<3-6 

<2-4 

MYSO W 33A (13) 


NH 2 CH 2 C 00 H '1 

< 0.3 

< 0.1 

MYSO W 33A (10) 



^ The environment f or which the upper limit was determined . 


*’ References: flVElsila et al. 199 

7. l2VSandford et al. 2001. 13 

1-Booeert et al. 2002a. 

f4V Vandenbussche et a 


199d. f5V< 

Sandford et al. 

1993 

. 161-Smith 

1991. 171-Smith et al. 

2011, (8)-Boudin et al. 

19981. 191-1 

I’ontoppidan et al. 

2008. 1101-C 

libb et al. 2004. IIII- 

van Broekhuizen et al. 

3 OO 5 I. fl2')-Schutte et al. ; 

999 

(13) 

-Bossa et al. 

2009 

” Detections toward comets (Mumma & Charnle 

V. 20111. 


Glycine 


Most species in Table [3] were studied because of their expected production in the 
gas phase or on grain surfaces. C 3 O 2 is a by-product of the energetic formation of 
CO 2 1 ^8.31) . Species for which, to the best of our knowledge, no detections were 
claimed or upper limits were determined include the electronically stable “gener- 
ation 0 ” spec i es wh ose formation was studied in the laboratory and are listed by 
Theule et al.1 (I2ni3l b HCN, CH 3 NH 2 , CH 2 NH, NH 2 OH, NH 2 CHO, CH 3 CHO, and 


HNCO. Other relevant species include O 3 , HOCN, and H 2 CO 3 . Energetic processing 
and thermal reactions efficiently produce more complex species, but these are increas¬ 
ingly hard to detect because of a lack of distinct absorption features, as evidenced by 
the high upper limit of glycine listed in Table [3l Millimeter wave spectroscopy of gas 
phase species is orders of magnitude more sensitive, and must be used to determine 
their abundances after sublimation off the grains. 


6.3 Abundance Variations 


A recurring conclusion in the literature is that the ice abundances relative to 
H 2 O are to the first order similar over a wide range of environments. As sample 
sizes and uncertainties have improved, clear trends of variations by factors of 2-10 
were discerned (Table [ 2 ]). Most of these relate to CO depletion and subsequent 
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Figure 8 : CO 2 abundances relative to H 2 O in a range of environments. For the 
comets, the abundance represents the ratio of gas phase p roduction rates (within 
the H 2 O sublimation radius at 2.7 AU: IOotsubo et al.ll 2012 lb For background stars, 


Taurus has distinctly lower abundances than other regions. 


chemical reactions 1 ^8.ip . A statistical study bv lOberg et al.l (j201lh revealed that 
A 1 h 2 o(CO) and -^H 2 o(C 02 ) are factors of ~3 and 1.5 larger toward LYSOs and 
the quiescent regions of isolated cores compared to MYSOs and Taurus background 
stars. Within samples of MYSOs and LYSOs, Ah 2 o(CH 30 H) and Ah 2 o(OCN“) 
(4.64 gm component) vary by an order of magnitude. CH 3 CH 2 OH or HCOOH fall 
in this category too. In quiescent regions, the observed Ah 2 o(CH 30 H) varies from 
<3% (Taurus) to ^12% in some isolat ed dense cores. 


The overview bv lOberg et al.l (j201lh does not include extragalactic and GO ices, 
which are particularly well studied in CO 2 . Figure [ 8 ] shows a comparison of Ah 2 o(C 02 ) 
in all environments with accurate measurements. It highlights the systematic differ¬ 
ences between environments, which are often rather small, but significant. MYSOs 
in the Galactic Disk have abundances that are typically a factor of 1.5-2 smaller com¬ 
pared to those in the LMC and compared to Galactic LYSOs. SMC MYSOs have low 
abundanc es when the HaO colu mn is low, but converge to the LMC values at higher 
columns ( Oliveira et ah . 2013 I l For background stars, the distribution is double- 
peaked. Low abund ances (~20%) are observed toward seven stars behind Taurus 


( Whittet et al.l . 120031 . while high abun dances (^40 %) are found for five lines of sight 


behind the S erpens dKnez et"^ . 200^1 and Lupus ( Boogert et ah . 2013) clouds, ari d 
the IC 5146 ( Whittet et ah . 2009 1. L429-C, and L483 cores ( Boogert et al.l . [2OI1I 1. 
The discriminator appears to be the line of sight extinction: all but one of the Taurus 
sightlines have Av< 25 mag, while all but one of the other sightlines have Av> 25 
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mag. It may thus be related to CO freeze out 1 ^8.11) . 

The 15.2 ^m CO 2 band profiles vary between envir onments as well. The 5 


component fits applied to LYSOs and Galactic MYSOs (iPontoppidan et al.l. 1200. 


were also applied to MYSOs in the GC CMZ ( An et al.l . 1201 il l and in the LMC 


( Seale et al.l . l201lh . The GC MYSOs have twice as much CH 30 H:C 02 and C0:C02 
ices relative to other environments (Fig. [9|). LMC MYSOs, on the other hand, have 
much stronger pure CO 2 components, indicating a larger degree of processing. Not 
shown here is the dominant C 02 :H 20 component, which is somewhat less abun¬ 
dant for the GC MYSOs and LMC MYSOs (~ 55%) compared to Galactic LYSOs 
(~ 70%). For the LMG sources, the sum of the pure CO 2 and G02:H20 ices is 
similar to the polar fraction toward LYSOs (jSeale et al.l . l201lh . suggesting that the 
pure CO 2 ices are “tapped” from the polar ices upon heating. For GG MYSOs, 
the CH 30 H:C 02 component makes up for the difference with LYSOs. Finally, some 
LYSOs have lower (G02:H20)/G02 fractions and higher Yh2o(C02 ), which is due to 
enhanced CO 2 formation in CO-rich environments (e.g., source SVS 4-10 in Fig. [5]). 

Variations of ice abundances relative to Vh are evident too. The CO depletion 


in ne arby quiescent clouds is incomplete (40-50%; IWhittet et al.lll985l : IChiar et al 
1994l i and the H 2 O abundance i s also a factor of ~2 below the values seen toward 
embedded YSOs ( Boogert et al. . 20131 ). At higher densities (10^ —10® 


cm ^) in cores 


PontoDPidan et al.ll2004l ) 


and envelopes all ice abundances increa se, but especially those o f CO and chemically 
related species (CH 3 OH, CO 2 , OCN” 

6.4 Isotopologues 


The ^^CO and ^^C02 isotopologues were securely detected toward interstellar 

/^^C02 ratios of 68-71 (two 
and 52-110 (thirteen sight- 


sightlines; 

Boogert et al. 

2002a: 

Pontoppidan et al. 

2003 

lines; 

Boogert et al. 

2000 

a) were derived, respectively. 


12 / 


C 02 /^®C 02 ratios in¬ 
crease with Galactocentric radius, in accordance with the elemental isotopic ratios. 
HDO was tentatively detec ted at an abundanc e of 2-22% relative to H 2 O toward 
LYSOs with edge-on disks ( Aikawa et al. . 2OI2I ). Otherwise, l ow upper limi t s wer e 


derived: <0.2-1% for in termediate mass an d massive YSOs ( Dartois et al. . 20031 ) 


and <0.5-2% for LYSOs ( Parise et ah . 20031 ). The relatively low de uterium fraction¬ 


ation of H 2 O is consistent with formation in an early cloud phase (jCeccarelli et al 

2OI4I : Ml . 
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Figure 9: The 15.2 ^um band of CO 2 shows distinct profile variations. Decompositions following Pontoppidan et al. ( 20081 ) show that the C 02 :CH 30 H 


and C02:CO components are particularly strong toward the GC MYSOs. LYSOs show a particular lack of C 02 :CH 30 H, while MYSOs in the LMC 
show an enhancement in pure CO 2 . Not shown is the C02:H20 component, for which median values are quoted in ^6.31 
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6.5 Interstellar Versus Cometary Ice Abundances 


Nearly all securely, likely, and possibly detected ice species outside the solar sys¬ 
tem are also detected toward comets (Table [2]) . Cometary ices are indirectly traced 
after out gassing and ions have recombined to their neutral counterparts. Interstellar 
NH^ should thus be compared to cometary NH 3 , OCN“ to HNCO, and HCOO“ to 
HCOOH. Overall, the median abundances of all C- and N-bearing species relative 
to H 2 O in comets are belo w those of LYSOs. For CH 4 , CH 3 OH, and CO, the dis- 
tribu tions overlap, howeve r (jOberg et al.l . 120111 1. Using the measurements in comets 
from Ootsubo et al. ( 2012 l 'l. a similar conclusion can be drawn for CO 2 ices (Fig. [ 8 ]): 
the median abundance is well below that of LYSOs (the upper quartile of comets 
is similar to the lower quartile of LYSOs), but there is some overlap. For CH 3 OH, 
it appears that the overlap is caused by the class of organic-enriched comets, which 
has members from both the Oo rt Cloud and the Kuiper Be lt, and has abundances 
similar to the median of LYSOs ( Mumma &: Charnlev . 201 li b There is no overlap in 
the distribution of NH 3 abundances, however. The discrepancy is even larger when 
the NH 3 and NH^ abundances toward LYSOs are added. 

Because cometary ices are traced indirectly using gas phase rotational transitions, 
more sensitive abundance measurements can be made compared to ice measurements 
toward LYSOs. This is evident in Table [3l which shows that upper limits for H 2 S, 
C 2 II 2 , C 2 He, and HNCO toward LYSOs are all above the detections for comets. 
Considering that HNCO in comets may directly trace OCN“ in LYSOs, the trend of 
low cometary C and N abundances are confirmed. The median is well below that of 
LYSOs, but there is some overlap. Finally, a number of species were detected toward 
comets, and no measurements are currently available for LYSO ices: HCN (0.08- 
0.5%), S 2 (0.001-0.25%), HOCH 2 CH 2 OH (0.25%), HCOOCH 3 (0.09%), HNC (0.003- 
0.05%), CH 3 CN (0.008-0.04%), HC 3 N (0.003-0.07%), and HaCS (0.05%). In view 
of the uncertain identification of a number of interstellar ice features, in particular 
those in the 5-8 pm wavelength region, a search for these ices is warranted. 


6.6 Elemental Budget 


The budget of chemical elements constrains the identification of detected features 
and helps direct searches for new volatile and refractory species. The securely-, 
likely-, and possibly-detected ices do not nearly account for the abundance of the 
heavy elements O, C, N, and S that is presumably available to be included in volatiles 
(i.e, excluding depletion in known refractory dust and gas phase volatiles). For O, 
~ 35% is “missing”, for C ~30%, for N 70-85%, and for S ~95%. 

The total amount of the available O included in H 2 O, CO, and CO 2 ices is similar 
in different q uiescent clouds (^26 % ). even though the abundances of the individual 


species vary ( Whittet et ah . 2009 I l Oberg et al. ( 2011 1 find a similar median value 


for MYSOs, and a somewhat larger median value of 34% for LYSOs, with a maximum 
of 61% in regions of high C O freeze out. T here is lingering uncertainty concerning 
putative additional carrier s (IWhittetl. I 2 OIC 1 L Sol i d O 2 could account for up to 10 % 
of O based on the limit in Vandenbussche et al. ( 1999 I l and the possibly identified 
species HCOOH for another few percent. 

Of the C budget presumably available for inclusion in volatiles, 27 and 14% is 
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taken up by known ices in LYSOs and MYSOs, respectively (jOberg et alJ . 120111 ). 
Maximum values of a factor of 2 larger are measured. Frozen PAH species are likely 
an important reservoir of (refractory) carbon in the ices, a t 18% of the volatile bud¬ 


get, a t least, using absolute band strength measurements (|Hardegree-Ullman et al. 


2014l i of neutral pyrene, but this ne eds to be re-a s sessed when the degree of ioniz ation 


of the frozen PAHs is constrained ( 

Keane et ah. 

2001 b: 

Bouwman et ah. 

2011 ). 

NH 3 , NH7, and OCN account for ~11% of the available N ( 

Oberg et ah. 

2011 ! 


It would be half of that if the identifications of NH 4 and OCN~ are incorrect. In 
some sightlines, up to 30% could be included in the ices. It is plausible that much of 
the “missing” N is in the form of solid N 2 , but this is extremely difficult to conhrm 
1 ^6.21) . The (uncertain) upper limits to the solid N 2 abundance (Table [3|) correspond 
to about 40% of the cosmic N budget. 

Finally, OCS and, if it s identification is conf irmed, SO 2 , contribute only 4% to the 
available sulfur budget (jPalumbo et al.l . Il997lb The upper limit to H 2 S (Table [3]) 
would account for another percent. 


7 GRAIN SIZE AND MANTLE THICKNESS 


Mantle growth is expected to proceed at a rate independent of the grain size, 
although the smallest grains (< 0 .002 /xm) will be i ce-les s as they become too hot 


upon stochastic photon heating ( Hollenbach et ah . 2009l h The elemental budget 


f ^ 6 . 6 p gives upper limits to the mantle thickness. For H 2 0-rich ices this is 0.005 gm 
assuming an MRN grain radius a distribution, at which point all O not included 
in gaseous CO and silicates is presumed to be incorporated into H 2 O. Subsequent 
complete freeze-out of CO increases the mantle thickness by a further factor ~ 2 . 
Such mantle growth can thus not significantly increase the sizes of the largest grains 
(a > 0.1 //m), but it can lead to a large fractional increase in the sizes of the abundant 
small grains that represent most of the surface area of the dust. Coag ulation of the 


sticky icy grains will lead to further growth in dense molecular clouds (jOrmel et al. 

2009( 1. 


The observed ice features place some constraints on the grain size and mantle 
thickness. There is strong evidence for grains that are small compared to the wave¬ 
length (a << \/2'k) with comparatively thick ice mantles. The peak optical depth of 
the H 2 O stretch mode is generally centered at 3.0 /xm (a << 0.5 ym), and the pro¬ 
files of the CO and CO 2 stretch modes are si ngle peaked. Thin mantled grains (core 
over mantle volume ratios larger than ~ 0.2; Tielens et ahl 1991 1 would have double- 
peaked profiles. In addition, spectrally resolved polarimetry of the 3.0 ym band 
shows a shift of the polarization peak to longer wavelengths i ndicative of dic hroic 


polarization (alignment of elongated grains by magnetic fields: IWhittetll201 la; An- 


dersson et ah, this volume) by grains smaller than ~0.1 /xm; (|Hagen et al 


19831 : 


Hough et ahl . 1996l l. To test the presence of ices on grains much smaller than this. 


observations of the electronic transition of H 2 O at 0.14 /xm would be needed. 

There is evidence for icy grains that are large compared to the wavelength as well, 
but with the excepti on of some clear-cut cases (e .g., abundant ~ 0.8 /xm icy grains in 
an LYSO disk; ^5.41 Terada &: Tokunaga 2012al l it has encountered problems. Mie 
scattering calculations show that grain sizes of up to ~0.5 /xm would be needed to 
explain the long wavelength wing of the 3.0 /xm ice band 1 ^4.11) . Such large grains are. 


































































Observations of the Icy Universe; ARAA 53 (2015; accepted), v. 05/05/2015 37 


however, inconsistent with the spectropolarimetric measurements mentioned above. 
Also, in models satisfying all obs ervational constrain ts (grain size distribution, ele¬ 
mental budget, ice band profiles) Smith et al. ( 19931 1 find that their observations of 
the 3.0 pm profile in Taurus can only be explained if ice mantles are preferentially 
present on specific grain populations (e.g., on larger grains or on silicate grains). 

Indirect evidence for a distinct population of much larger, micron-sized icy grains 
comes from diffuse s cattered light observed toward dense cloud cores (“coreshine”; 


Andersen et al.ll2014l '). Such grains were also suggested as an al ternative exp lanation 


of the long-wavelength wing of the 4.67 pm CO ice band ( 114.21 lDartoisll2006l). Larg e 
grains have been suggested as a source of the “missing” oxygen ( ^6.61 Jenkins 2009lb 
because for radii larger than a few micron, the infrared absorption profiles become 
too wide and shallow to be detected. In light of the small variations of the 3.0 pm 
band profile ( 114.11) . these suggestions need further study, including an assessment of 
the effects of grain coagulation on the band profiles. 


8 CONSTRAINTS ON ICE EVOLUTION MODELS 

The commonalities and differences in the ice abundances and band profiles observed 
over a wide range of environments provide strong constraints on models of ice for¬ 
mation and evolution. We do not review the models here, but rather summarize the 
observational constraints. Figure [TOl illustrates the probable scenario. 

8.1 Accretion and Grain Surface Chemistry 


The observed ice abundances are generally in good agreement with models of 
accretion from the gas phase followed by grain surface chemistry. In particular, 
the efficiency of h ydrogenation as a resu lt of the mobility and tunneling capacity 
of H atoms (e.g., Tielens fc Hagen 19821') is confirmed by the observations. The 
dominance of H 2 O ice and mixtures with H 2 O in all Galactic and extragalactic 
environments points to a phase of molecule formation on cold grain surfaces at 
which the accreting gas is atomic H- and 0-rich. This is a phase of relatively low 
densities (n > 10^ cm~^) early in the cloud evolution (H/H 2 ratios decrease rapidly 
with density; Hollenbach et all 1971 1. The ice formation threshold of Av=1.6 mag 
near local cloud edges is understood as the onset of rapid ice mantle growth after 
a monolayer has formed and the effect of photodesor ption decreases with inc reasing 
extinction by dust into the cloud (Ayoc ln{Go/n); Hollenbach et al. 2 OO 9 I 1 . This 
phase also includes the formation of CH 4 , CO 2 (polar component, likely via CO-I-OH; 


loDDolo et al.ll201ll l. and NH 3 which are all intimately mixed with H 2 O. Variations 


in CO 2 /H 2 O ratios in this early phase reflect variations of the CO accretion rates. 


more 0 can be included in H 9 O (Whittet et ah. 2009: 

Oberg et ah 

201 Ih It mav 

explain the possiblv C09-free region in the SMC (Oliveira et ah. ^ 

1013') and CH/i- 

rich region in Galactic cloud edges (Oberget ah. 2011 

). On the other hand, CO 

accretion rates increase at lower temperatures and at higher densities, leading to 


enhanced formation of CO 2 (up to a factor of 2 ), still well mixed with the co-formed 
H 2 O species. This phase is CHsOH-poor and has a relatively low CO freeze out 
(< 50%). At the same time, the presence of the 6.85 ^m feature and some absorption 

























































38 


Boogert, Gerakines, Sz Whittet 


underlying the 6.0 fim band (Fig. [ 6 |) show that these early ices are already quite 
com plex. They are formed on a short time scale (10®/n r\j 10 ® yr at n = 10 ^ cm 
e.g., Hollenbach et ahlboogl l compared to the cloud life time (~ 10 ^ yr) and compared 
to the formation time of new species by interstellar CRs or CR-induced UV fields 

(H). 

Strong enhancements (up to a factor of 10) of the abundances of CH 3 OH, OCN“, 
the carriers of the 7.24 and 7.41 pm features (CH 3 CH 2 OH or HCOOH and CH 3 CHO 
or HCOO“), and the apolar components of the CO and CO 2 ices highlight the 
subsequent phase of ice formation. It is the episode of “catastrophic” CO freeze 
out at low temperatures and high densities (n > 10 ® cm“^), when the gas phase 
H/CO ratios are drastically enhanced, even though the H/H 2 ratios are reduced. 
This then leads to H 2 CO and CH 3 OH on a short time scale of ^ few x 10^ years, as 


demonstrated in Monte Carlo simulations (ICuppen et al 


of other CO-derived species, including HNCO (jFedoseev et al 


2 OO 9 II . and l ikely to a host 
2 OI 5 I 1 and CO 2 (i 


m 


the apolar component). The best observational evidence of this phase comes from 
the measurements of ice abundance gradient s in dense cores and Class 0 LYSOs 
( Pontoppidan et ah . 2004; Pontoppidan . 20061 '). The quiescent regions in the Taurus 
cloud have not entered this phase yet (Fig. [S]), in contrast to a subset of dense cores, 
prestellar cores, and LYSOs (e.g., S VS 4-10 in Fig. [5] has a C O 2 /H 2 O ratio of 50% 

Clearly, a complex 


and a strong CO:C 02 component; iPontoppidan et al. 


iooib 


mixture of ices is formed before stars form. Apparently, the GC MYSOs (prominent 
apolar CO 2 and CH 30 H:C 02 components) have gone through this phase, but the 
LMC MYSOs have not (weak apolar CO 2 and CH 30 H:C 02 components and low 
solid CO abundances). We speculate that high and low metallicities in these distinct 
environments, respectively, play a role. High metallicities reduce the H/CO ratios 
and increase extinction and thus CO freeze out. 






































1. Polar Phase: gas H/CO large 
T>20 K, n>103 cm-3, Av>1.5 


Grain Surface Chemistry 

2. Apolar Phase: gas H/CO small 
T<20 K, n>104 cm-3, Av>3 


3. CH 3 OH Phase: gas H/CO large 
T<20 K, n>105 cm-3, Av>9 


CR/UV ^ 


h| 
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Thermal Processing 


4 . Distillation Phase: 

T>20 K 
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T~30-77 K 
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T>90 K, high UV ,jv 

CR 



CO 

CO 


Figure 10: Schematic view of the observationally constrained composition of icy grain mantles, highlighting the physical and chemical processes at 
play. Six phases are shown, which could represent progressive evolutionary phases in quiescent clouds (phases 1-3) and YSO envelopes and disks 
(phases 3-6). Gas phase species participating in grain surface reactions are indicated in red. In phases 1-5, the ices experience CR and secondary 
UV fields that only minimally affect the main ice species. As the ice temperature rises (phases 4 and 5), thermal reactions or reactions involving 
radicals may proceed. Stage 6 is least observationally constrained. The relative grain and mantle sizes in the first 5 stages are realistic (maximum 
mantle thickness of 50 A and minimal grain core sizes of 20 A). After the initial ice layers are formed (phase 1), grain coagulation accelerates. 
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8.2 Thermal Processing 


Observations show that after the formation of the bulk of the ices through accretion 
and grain surface chemistry, the evolution of the ices in the circumstellar envelopes 
and disks of YSOs is dominated by thermal processes (S|). In summary, sublimation 
of apolar ices occurs at temperatures of ~ 20 K. As CO is “distilled” out of the CO 2 
ices, this develops the 15.11 and 15.26 gm absorption peaks of pure CO 2 . Segregation 
of CO 2 from H20-rich ices occur s at higher temperat ures, starting at 30 K for the 


surface l ayers and thin mantle s ( Oberg et ah . 2009bl k and up to ~77 K for thick 


mantles (|Boogert et al.l . l2000al ) at which the H 2 O crystallization process starts as 
well. This is observed by a long-wavelength shift and narrowing of the 3.0 /rm band, 
and again by the development of the 15.11 and 15.26 gm absorption peaks. The 
“pure” CO 2 will sublimate at temperatures in the range of ~45 K for “distilled ices”, 
up to the H 2 O sublim ation temperature o f ~90 K for segregated ices, depending on 
the mantle thickness (IFavolle et al. 1 . 120111 1. Sublimation of t he polar, H20-r i ch ice s 
is traced by increased gas/solid state column density ratios (jBoonman et al.l . l2003lb 
The strength of the C4 component of the 6.85 /rm featu re is apparently enha nced at 


temperatures above the H 2 O sublimation temperature ( Boogert et ah . 2008l b 


The sublimation, segregation, and crystallization temperatures quoted in this re¬ 
view apply to the long time scales in the interstellar environments. At any temper¬ 
ature, the probability that the energy barrier for these processes is surmounted is 
set by the Boltzmann d i stribut ion. This translates to time scales of i/~^e^baTTieT/T 


(jTielens fc: Allamandolal . Il987bl b For example, for H 2 O, with an energy barrier of 
.^'barrier = 5500 K and a frequency = 5 x 10“^^ s“^ (0-H bend vibration mode), 
an observed sublimation temperature of 180 K at a laboratory time scale of minutes 
corresponds to 90 K at space time scales of 10,000 years. These arguments were used 
to link the range of pure CO 2 peak strengths in a sample of MYSOs to a time scale 



ice sublimation radius of 1000 AU, a crystallization and segregation region between 
100 and 130 AU, and polar ice sublimation within 100 AU. In circumstellar disks 
surrounding similar stars, the observed ice crystallization 1 115.4P is also thought to 
happen at radii > 1 00 AU, but only in their super-heated optically thin surfaces 
( Chiang et al.l . 200 ih . 

For objects of very low lu minosity, e n velop e heating is insufficient to cause crys¬ 
tallization and segregation (IKim et al.l . [2013). The double peaked CO 2 profile is 
detected toward <0.1 Lq YSOs, however. It is suggested that the apolar CO 2 com¬ 
ponent is responsible for that, in combination with heating spikes caused by episodic 
accretion events. During long intervals between the accretion events, CO 2 is formed 
in the CO-rich ice by reactions with OH. The heating event then sublimates the CO, 
and leaves pure CO 2 behind. In a separate study, a CO 2 bending mo de consisting 


of ne arly pure CO 2 is found toward the 1 Lq YSO HOPS 68 (Fig. O jPoteet et al 


2013lb This may also be interpreted as a result of episodic heating, but in this case 
the event is much stronger (up to 100 Lq, as in FU Ori objects) and all ices are 
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sublimated and then recondensed in the order of their condensation temperature. 
Fortuitous timing and line of sight geometry inside a flattened envelope are required 
to trace this particular ice component. 

The common process of ice heating most likely affects the chemistry as well, al¬ 
though from an observational point of view there is little direct evidence for this. A 
candidate feature produced upon heating is the “C5” component (Fig. [6]). We spec¬ 
ulate that it might be formed when the radicals, created by energetic processes over 
the cl oud and envelope life time, become more m obile and react at higher tempera¬ 
tures ( Garrod et ah . 2008 : Vasvunina et ah . 2014l b Finally, purely the rmal reactions , 


without the help of energetic processe s , may well be quite common (jSchutte et al. 


witnout tne neip or energetic processe s , may ' 
I 993 I : Bossa et ah . 2009 : Theule et al. . 2013). 


The formation of NH^ and OCN 


(proposed) carriers of the 6.85 and 4.62 pm features, proceeds in f act by an acid- 


base reaction of NH.s with HNCO a t temperatures as low as 15 K (jRaunier et al. 


2004bl : Ivan Broekhuizen et al.l . [2004l i . This is consistent with NH^ being the carrier 
of the 6.85 pm feature, as it is observed in quiescent sightlines. HNCO cannot be 
the only acid involved in this reaction, however, as the OCN“ abundance is an order 
of magnitude below that of NH4 . 


8.3 Energetic Processing 


Historically, it has long been held that the energetic processing of ices in space plays 
a key role in the production of new ice species and specifically in th e transformation 
from simple to more complex components (e.g., iGreenberd Il973l ). Ices in dense 
clouds experience C R fluxes that are represented by 1-MeV protons with a flux of 

and references therein). The interstellar UV 


1 proton cm ^ s ^ 


(iMoore et al 


2001 


field of ~ 10® Go eV cm ® s ^ (with Go = 1 in the local ISM) rap idly declines with 


depth into the cloud, regulating the on-set of ice formation (e.g., iHollenbach et al 


2009l b Within clouds, the fluxes of UV photons are much lower, generated by GR 
destruction of H 2 and corresponding to ~ 10^ eV cm“^ s“^. In regions of star 
formation these fluxes may be strongly enhanced in both space and in time. 

A vast literature of laboratory studies of energetic pr ocessing exists, which has 
been particularly successful in solar system studies (e.g., Paranicas et al. 2009l b In 
contrast to grain surface chemistry and thermal processing, the significance of en¬ 
ergetic processing of ices outside the solar system has not yet been unambiguously 
demonstrated, however. Observational tests of energetic processing must address the 
similarities and differences between laboratory and observed interstellar abundances, 
as well as the spectral profiles of the produced species of interest and all subsequent 
chemical (by-)products. 

Many ice absorption features have been studied in the context of energetic pro¬ 
cessing, but most attention was paid to those of GO 2 , as well as those at 4.62 and 
6.85 pm. GO 2 is formed readily by the energetic processing of most carbon-bearing 
molecules in the prese n ce of H 2 O, including pure H 2 O ice mantles on carbon grains 
( Mennella et ah . 2004l f . IoddoIo et al. (2013) have explored this to the full extent by 
fitting the 15.2 pm GO 2 bending mode toward MYSOs, LYSOs, and dense clouds 
using only CO 2 produced by processing with Galactic GRs. Gomparing these results 
to the empirical decomposition (Fig. [5|), the broad polar component is fitted with 
GO 2 produced from irradiated H 2 O on carbon grains, the shoulder from irradiated 
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CH 3 OH, the apolar components from irradiated CO and CO:N 2 ices, and the pure 
CO 2 component by heating of these ices. At the observed abundances of CO, H 2 O, 
and CH 3 OH, the CO 2 profiles and abundances are reproduced after ~ 10^ yr for 
most targets. Similarly long time scales were found fo r the production of CO 2 by 


the CR-generated UV field (see above) in dense clouds (jWhittet et al.l . Il998li . Thus, 


overall, after the initial rapid formation of CO 2 and H 2 O by grain-surface reactions 
(in about 10 ^ yr; see 118.11) . additional CO 2 may be formed by particle bombard¬ 
ment during the millions of years in YSO environments, or on shorter time scales 
for regions exposed to stronger UV and particle fields. Not all implications of this 
model have been fully investigated, however. Most important ly, the effect of a suite 


of by-products, in particular carbon-chain oxides (C 3 O, C-^Oo: iGerakines et al. 


1996 


Gerakines fc: Moord l2001al : iPalumbo et al.l I 2 OO 8 I I on other spectral regions has not 
yet been fully assessed. Also, no distinct enhancements of GO 2 abundances were 
observed toward YSOs (Fig. [8]), although the situation is complicated by the suscep¬ 
tibility of GO 2 itself to destruction by irradiation. 

The 4.62 gm XGN feature was seen for a long time as an indicator of ener¬ 
getic processing, as laboratory UV photolysis and ion bombardment experiments 
on mixed GO and NH 3 ices successfully reproduced this absorption feature with 
OGN“ ( 116.11). Similar experiments c laimed an identification of the 6.85 Atm featur e 


with NH^ ( Schutte fc Khanna . 2003l l and H 2 NGONH 2 ( 


urea; 


Raunier et al 


>nn4alb 


Both features can also be reproduced through purely thermal acid-base reactions in 
ices, however and in many sightlines this is the more likely process 


8.4 Connection with Cometary Ices 


The survival of the protostellar envelope ices upon accretion into the protoplane¬ 
tary disk is an open question. Depending on the radius to the central s tar, full or par- _ 

tial d estruction of ices and subsequent reformation may have occurred (|Mumma &: Gharnlev 
201 ll . and references therein). The Stardust mission revealed that refractory dust 


may have formed close to the early Sun, moving radially outward, accumulating ices 
along the way and thus dilu ting any presolar components in the outer, comet forma¬ 
tion zones ([Brownlee! . l2014l b Ice observations provide constraints to this discussion, 
but have not been conclusive. There are distinct differences between the composition 
of cometary and LYSO envelope ices 1 ^6.51) . G- and N-bearing species in comets are 
under-abundant on average. For G, the abundance distributions do overlap, which 
may imply that the early solar system resembled an LYSO with less G in the ices 


(jOberg et al.l . [2011]). Low cometary GO and GO 2 may also reflect sublimation of 
the most volatile, apolar component s upon accretion from the envelope to the pro¬ 


toplanetary disk ([Visser et al 


J ) onent s upon accretion irom tne envelope to tne pro- 
, 2 OO 9 I L Neither argument holds for N, however: the 


NH 3 abundance distributions do not overlap and NH 3 is only present in the polar 
(least volatile) ices. Direct observations of ices in the disk mid-plane, where comets 
are formed, are not possible. Ices in the disk surface may be linked to the mid¬ 
plane by turbulence, however. High spatial resolution observations have shown that 
these ices experience thermal processing 1 ^5.41) . Further observations are needed to 
search for evidence of energetic processing, but chemical models including vertical 
disk turbulence have indicated that the ice composition through out the disk may be 
significantly affected by energetic processes in the disk surface ([Giesla &: Sandford . 
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20121 ') ■ Radial transport might lead to such effects as well. 


9 SUMMARY POINTS 


Infrared observations have revealed the icy Universe (Fig. I10|) : 

1 . H 20 -based ices trace dense clouds (n > 10 ^ cm“^) at depths corresponding to 
Ay> 1.6 and temperatures below 90 K. These are mixed ices, containing at 
least CO 2 , NH 3 , CH 4 , and the carrier of the 6.85 pm feature (NH 4 ?). 

2 . CO freeze-out occurs deeper in dense clouds 3), forming an apolar ice 

layer, and it does so “catastrophically” at 9 and n > 10 ^ cm“^, initiating 
a CO-based chemistry on the ices forming a third, CHsOH-rich, layer including 
also H 2 CO, newly formed CO 2 , and possibly HNCO. 

3. Heating of the ices is commonly observed and the spectroscopic signatures 
provide diagnostics of the thermal history of dense environments. Crystalliza¬ 
tion signatures in the 3.0, 45, and 63 pm H 2 O bands, and the position of the 
6.85 pm feature trace temperatures > 77 K. The double-peaked 15.2 ^m CO 2 
band traces a combination of such high temperatures (“segregation”) and tem¬ 
peratures >20 K as a result of CO sublimation (“distillation”). The ratio of 
the apolar/polar components of the 4.67 pm CO band itself traces the latter 
conditions as well. 

4. Most observed ice features are consistent with molecule formation on cold 
grain surfaces. The observed ice heating may well induce chemical reactions, 
possibly involving radicals created by energetic processes over long time scales. 
The observational evidence for this, possibly features in the 5-7 pm region, is 
weak, however. 

5. Variations of the CO 2 /H 2 O column density ratios and CO 2 ice composition 
in local versus GC and extragalactic environments likely reflect metallicity 
differences. 

6 . Median abundances of C and N-bearing cometary ices are below those of the 
envelopes of LYSOs. The abundance distributions do overlap, however, and 
inheritance of cometary ices from the ISM and envelope phases can thus not 
be excluded. 

7. Laboratory spectroscopy has been crucial in the identification and in the de¬ 
termination of the chemical and physical history of the observed ice features. 

10 FUTURE ISSUES 


1. Recent high spatial resolution observations of circumstellar disks have shown 
that the ice composition and processing history can be traced on scales relevant 
for p l anet and comet formation, both spectroscopically (e.g.. ISchegerer fc Wolf 


2016 : 


20091 : 


Terada &: Tokunagall2012a ') and using scattered light imaging (e.g., Honda et al. 


Debes et al.ll2008l L Future high angular resolution ground (Thirty Meter 


Telescope, Extremely Large Telescope) and space-based (James Webb Space 
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Telescope; JWST) infrared facilities thus promise great advances of ices studies 
at these scales. 

2. The same new facilities will enable ice studies in extreme environments (strong 
radiation fields, high or low metallicities), such as Galactic and extragalactic 
MYSOs at larger distances and in crowded fields. Observed feature profile 
and depth variations may further constrain identifications as well as molecule 
formation pathways, including the much sought after evidence for energetic 
processing. 

3. The increased sensitivity of new facilities will drastically enhance studies of 
quiescent clouds at very low and high extinctions, which are needed to study 
the ice formation process and its interplay with the grain growth process. 

4. The continuous 1-29 fim wavelength coverage provided by JWST at spectral 
resolving powers R of up to 3000 will provide strong constraints on the iden¬ 
tification of new ice species. Supporting laboratory experiments are essential 
for this. 

5. New instruments on platforms such as the Stratospheric Observatory for In¬ 
frared Astronomy (SOFIA) and the Space Infrared Telescope for Cosmology 
and Astrophysics (SPICA) are needed to give access to the ice lattice modes 
(> 25 /um), which are particularly powerful tracers of the ice composition and 
thermal history. As they appear in emission, the ices may be mapped. 

6. The good understanding of the main ice bands allows them to be used as tools 
in the characterization of new sources that will be found with new infrared 
facilities. 

11 RELATED RESOURCES 

Infrared transmission spectra and optical constants of astrophysically relevant pure 
ices and mixtures may be found at the following web sites; 

• Astrophysical and Astrochemistry Laboratory at NASA Ames Research Cen¬ 
ter: 

http://astrochem.org/db.php 

• Experimental Astrophysics Laboratory at Catania Astrophysical Observatory: 
http://www.ct.astro.it/lasp/ 

• Sadder Laboratory for Astrophysics (Leiden): 
http://www.laboratory-astrophysics.eu/ 

• Cosmic Ice Laboratory at NASA Goddard Space Flight Center: 
http://science.gsf c.nasa.gov/691/cosmicice/ 

• Grenoble Astrophysics and Planetology Solid Spectroscopy and Thermody¬ 
namics database (GhoSST): 

http://ghosst.osug.fr 

• Jena - St.Petersburg Database of Optical Constants (mostly references to lab¬ 
oratory spectra): 

http://www.astro.uni-j ena.de/Laboratory/Database/jpdoc 

• Virtual Atomic and Molecular Data Centre (VAMDC); 
http://portal.vamdc.org 
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Computer programs for calculating the effect of grain shape and size on ice ab¬ 
sorption band profiles can be found in Bohren fc Huffman ( 198ill ) and are also made 
available by B. T. Draine at 


http://www.astro.princeton.edu/~draine/scattering.html 
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